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ABSTRACT 


After  a  brief  review  of  the  history  of  material  failure  studies  and 
of  measurements  of  stored  energy  in  plastically  deformed  bodies,  energy  and 
stress  and  strain  relations  are  developed  for  both  elastic  and  plastic 
deformation  in  uniaxial  strain.  The  physical  effects  of  plastic  deformation 
are  discussed  and  foundations  are  laid  for  thermodynamic  calculations  of 
plastically  deformed  materials.  Calculations  of  both  thermodynamic  and 
mechanical  parameters  are  described  for  uniaxial  strain  with  and  without 
energy  storage  in  internal  strains.  Conditions  of  uniaxial  strain  are  found 
to  yield  new  thermodynamic  relations  which  are  useful  for  synthesizing  con¬ 
stitutive  relations. 
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I .  INTRODUCTION 


From  the  earliest  times  men  have  been  concerned  with  the  conditions  of 
failure  of  structural  materials.  Little  is  known  of  early  Egyptian,  Greek  and 
Roman  knowledge  of  the  subject,  though  it  is  evident  from  their  extensive  con¬ 
structions  that  they  had  some  effective  working  rules.  During  the  Renaissance 
there  began  to  be  recorded  systematic  attempts  to  study  and  develop  the  con¬ 
cepts  of  material  failure.  Leonardo  da  Vinci  describes  an  apparatus  for  testing 
the  strength  of  iron  wire.  In  a  note  on  "Testing  the  Strength  of  Iron  Wires  of 
Various  Lengths"  he  says,  "...The  object  of  this  test  is  to  find  the  load ’an  iron 
wire  can  carry.  Attach  an  iron  wire  two  braccia  long  to  something  that  will 
firmly  support  it,  then  attach  a  basket  or  any  similar  container  to  the  wire  and 
feed  into  the  basket  some  fine  sand  through  a  small  hole  placed  at  the  end  of  a 
hopper.  A  spring  is  fixed  so  it  will  close  the  hole  as  soon  as  the  wire  breaks. 
The  basket  is  not  upset  while  falling,  since  it  falls  through  a  very  short  dis¬ 
tance.  The  weight  of  sand  and  the  location  of  the  fracture  of  the  wire  are  to 
be  recorded.  The  test  is  repeated  several  times  to  check  the  results.  Then  a 
wire  of  one-half  the  previous  length  is  tested  and  the  additional  weight  it 
carries  is  recorded;  then  a  wire  of  one- fourth  length  is  tested  and  so  forth, 
noting  each  time  the  ultimate  strength  and  the  location  of  the  fracture." 
(Timoshenko,  1953)  It  is  particularly  interesting  to  note  that  even  here,  in 
the  fifteenth  century,  the  concepts  of  graduated  tests,  systematic  recording  of 
observations,  reproducibility,  and  variation  of  a  single  parameter  while  holding 
others  constant  were  fully  developed. 

Many  of  the  most  famous  names  in  science  were  involved  with  problems  of 
material  failure:  Galileo,  Euler,  Lagrange,  Coulomb,  Poisson  and  Kelvin,  to  name 
but  a  few.  The  basic  concepts  of  plasticity  were  developed  by  Barre  Saint-Venant 
and  described  in  several  papers  during  the  latter  half  of  the  nineteenth  century. 
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This  was  a  period  of  intensive  activity  in  mechanics  and  one  subject  of  interest 
was  the  search  for  a  suitable  criterion  for  failure  of  ductile  materials,  dames 
Clerk  Maxwell  described  the  distortion  energy  concept  of  failure  in  a  letter  to 
Lord  Kelvin  in  1856.  Although  criteria  based  on  this  concept  are  widely 
accepted  today,  it  did  not  appear  in  print  until  the  early  part  of  this  century, 
when  it  was  developed  independently  by  M.  T.  Huber  and  R.  von  Mises. 

Since  1934  when  E.  Orowan,  M.  Polanyi  and  G.  I.  Taylor  independently 
developed  the  concept,  and  especially  since  1950,  failure  theories  of  ductile 
solids  have  been  based  primarily  on  dislocation  behavior.  A  great  body  of 
information  on  dislocations  and  their  properties  and  behavior  has  grown  up 
during  the  last  twenty  years,  and  it  is  realized  that  the  continuum  concepts  of 
failure  are,  in  a  microscopic  sense,  ill  founded.  It  is  not  possible  to  dis¬ 
place  them,  however,  because  the  theory  of  dislocations  is  incomplete,  because 
it  is  much  more  complicated  than  continuum  theories,  because  continuum  theories 
work  quite  well,  and  because  there  have  been  inadequate  and  only  partially 
effective  efforts  to  bridge  the  chasm  which  separates  the  microscopic  behavior 
of  dislocations  and  the  macroscopic  behavior  of  real  materials. 

Because  of  this  gap  between  microscopic  understanding  and  macroscopic 
behavior,  there  has  been  an  intensive  parallel  effort  to  develop  and  refine 
continuum  concepts  of  plasticity  and  yield  during  recent  years. 

Until  the  1940's  there  was  little  concern  about  the  thermodynamics  of 
plastic  strain.  In  1948  Carl  Eckart  published  an  important  paper  on  the  thermo¬ 
dynamics  of  anelastic  strain,  and  Tolman  and  Fine  discussed  applications  of  the 
Onsager  theory  of  irreversible  thermodynamics  to  problems  of  plasticity.  These 
were  followed  by  Bridgman's  paper  of  1950,  an  article  by  Reiner  in  1958,  and 
then  by  a  spate  of  articles  in  the  1960's,  of  which  perhaps  the  best  known  are 
those  by  B.  Coleman  and  his  coworkers,  Meixner,  and  by  G.  A.  Kluitenberg. 
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Except  for  Bridgman's  paper,  most  of  these  have  been  concerned  with  general  and 
quite  abstract  formulations.  Indeed  it  is  quite  often  very  difficult  to  ferret 
out  the  underlying  physical  assumptions  on  which  the  work  is  based. 

There  appear,  then,  to  be  two  needs  for  current  work.  One  is  to  formulate 
explicitly  and  sharply  physical  foundations  for  a  thermodynamic  theory  of 
plasticity  which  can  be  grasped  by  other  than  trained  mathematicians  and  then  to 
apply  these  to  a  relatively  simple  yet  experimentally  realistic  physical  situation. 
The  other  is  to  join  the  continuum  thermodynamic-mechanics  developments  to 
dislocation  theory  in  a  manner  which  may  be  useful  to  both  subjects. 

In  the  present  report  some  progress  is  made  toward  each  of  these  goals. 
Further,  since  the  question  of  stored  energy  in  plastic  work  is  everpresent  and 
basic  to  the  subject,  a  brief  review  of  the  state  of  stored  energy  measurements 
is  given. 

The  principal  results  of  the  theoretical  development  are  that 

(1)  a  point  of  view  is  developed  that  assumes  the  thermodynamics  of  plastic 
deformation  to  be  essentially  reversible  if  the  correct  variables  are  chosen; 

(2)  theoretical  calculations  and  experiments  in  one-dimensional  strain  introduce 
density  as  a  variable,  and  this  leads  to  phenomena  and  thermodynamic  relations 
whose  utility  goes  beyond  that  found  for  conditions  of  uniaxial  stress; 

(3)  a  combination  of  thermodynamics  and  dislocation  theory  is  fruitful  for  both 
fields,  but  the  required  relations  are  not  completely  developed  here. 
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II.  EXPERIMENTAL  MEASUREMENTS  OF  STORED  ENERGY  OF  DEFORMATION 


Methods  for  measuring  the  stored  energy  in  a  plastically  deformed  solid 
may  be  put  into  two  categories: 

(1)  single-step  methods  in  which  all  the  necessary  measurements  are  made  during 
the  deformation; 

(2)  two-step  methods  in  which  the  stored  energy  is  measured  following  deformation. 
2.1  Single-Step  Methods 

These  methods  are  based  on  a  direct  application  of  the  first  law  of 
thermodynamics.  The  stored  energy  E$  is  the  algebraic  sum  of  the  work,  W, 
and  the  heat,  Q.  Work  may  be  evaluated  from  stress  and  strain  measurements.  The 
heat  effect  is  determined  either  from  temperature  rise  of  the  specimen  deformed 
under  adiabatic  conditions  or  by  calorimetry. 

The  first  successful  attempt  to  measure  the  stored  energy  in  a  material 
undergoing  a  deformation  was  made  by  Hort  (1906,  1907).  In  his  experiments, 

Hort  conducted  tensile  tests  in  iron  in  a  water  calorimeter.  The  calorimeter 
was  calibrated  by  electric  heating  and  thermoelastic  cooling.  For  true  strains 
ranging  from  0.034  -  0.109,  the  measurement  of  stored  energy  in  iron  gave  values 
of  2.2  to  5.1  cal/g  atom. 

Farren  and  Taylor  (1925)  determined  the  heat  evolution  associated  with 
nearly  adiabatic  extension  of  polycrystall ine  specimens  of  steel,  copper,  and 
aluminum,  and  of  a  single  crystal  of  aluminum.  The  samples  were  pulled  for  a 
short  time  and  work  was  obtained  from  a  recorded  stress-strain  curve.  The  heat 
evolved  was  computed  from  temperature  rise  in  the  specimen  measured  with  a  thermo¬ 
couple  placed  inside  an  axial  hole  in  the  specimen.  In  1934,  Taylor  and  Quinney 
repeated  these  measurements  by  deforming  steel  and  copper  under  torsion  to  obtain 
higher  strains.  Mechanical  work  was  evaluated  as  in  their  previous  work  while 
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heat  was  measured  both  by  a  thermocouple  and  by  a  water  calorimeter.  The  two 
methods  gave  similar  results. 

Maier  and  Anderson  (1934)  were  the  first  to  measure  energy  stored  during 
wire  drawing  by  the  single  step  method.  Their  working  metals  were  aluminum  and 
copper. 

The  measurement  of  stored  energy  under  compression  was  initiated  by 
Kunin  and  Senilov  (1936).  They  compressed  copper  by  stages  in  a  testing  machine; 
work  was  measured  by  special  "crusher"  gauges.  Epifanov  and  Rebinder  (1949) 
determined  the  energy  stored  in  chips  of  aluminum  formed  by  drilling  specimens 
inside  a  calorimeter.  The  difference  between  the  energy  input  to  the  drill  and 
heat  evolved  was  considered  to  be  the  energy  stored  in  the  chips. 

Eugene  (1953,  1954)  determined  stored  energies  in  lead,  copper,  and 
steel  by  using  a  freely  falling  weight  to  deform  the  specimens.  He  obtained 
stored  energy  from  the  difference  between  the  known  energy  of  the  weight  striking 
the  specimen  and  the  heat  effects  as  measured  by  its  temperature  rise.  He  was 
able  to  find  very  small  amounts  of  stored  energy  i.e.  less  than  1%  compared  to 
the  average  of  around  10%  reported  by  earlier  workers. 

Williams  (1960)  developed  a  device  in  which  specimens  were  deformed  by 

freely  falling  hammers  of  tungsten  carbide.  Since  the  period  of  contact  is  less 
-3 

than  2  x  10  sec,  the  deformation  induced  in  the  sample  is  assumed  to  be  adiabatic. 
The  sample  is  supported  by  a  single  thermocouple  which  measures  the  temperature 
rise.  The  work  term  is  difficult  to  evaluate  because  of  the  mechanical  losses 
due  to  suspension  of  the  hammers.  Heat  transfer  between  the  hammers  and  the 
sample  is  small  because  of  the  short  interval  of  contact  between  them.  Heat 
loss  due  to  interaction  with  the  surroundings  is  minimized  by  performing  the 
experiments  in  vacuum.  Williams  (1961,  1962,  1963)  has  determined  the  stored 
energy  in  copper  and  copper-silver  alloy  by  this  method. 
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All  of  the  measurements  described  above  were  carried  out  at  room  tem¬ 
perature. 

Erdmann  and  Jahoda  (1964)  designed  a  device  for  deforming  metal  wires 
at  low  temperatures  and  for  simultaneously  measuring  thermal  conductivity, 
electrical  conductivity,  and  heat  liberated  by  the  deformation.  The  experiment 
is  beset  with  the  usual  difficulties  associated  with  low  temperature  experi¬ 
ments  which  leads  to  some  uncertainty  in  data  reduction. 

Williams  (1963)  also  developed  yet  another  apparatus  in  which  the  heat 
developed  during  the  deformation  of  a  sample  under  tension  is  measured  by  the 
amount  of  vaporization  of  a  working  fluid.  The  method  is  generally  applicable 
from  room  temperature  to  77°K.  Wolfenden  and  Appleton  (1963)  have  developed 
a  method  similar  to  the  calorimetry  of  Williams.  The  fluid  used  is  liquid 
nitrogen.  These  authors  determined  the  stored  energy  in  both  single  and  poly¬ 
crystals  of  copper  and  aluminum. 

0.  W.  Dillon  (1962,  1966,  1968)  in  a  series  of  papers  determined  the  heat 
generated  during  torsional  oscillations  of  aluminum  and  copper  tubes  and  the 
heat  generated  during  the  propagation  of  plastic  waves  in  aluminum.  This 
apparatus  is  based  on  the  design  of  Farren  and  Taylor  (1925). 

2.2  Two-Step  Methods 

In  these  methods  a  material  is  deformed  plastically  by  any  suitable 
process.  The  energy  spent  in  deformation  may  or  may  not  be  measured.  The 
stored  energy  is  determined  in  a  subsequent  step  by  measuring  the  difference  in 
the  thermal  behavior  of  the  deformed  and  a  reference  specimen  in  attaining 
identical  final  states.  The  difference  in  these  heat  effects  is  the  stored 
energy  of  cold  work. 
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Two  step  methods  may  be  divided  into  two  general  classes. 

(1)  Annealing  methods. 

(2)  Reaction  methods. 

(1)  Annealing  methods  may  be  further  subdivided  into: 

(a)  Isothermal  annealing,  where  the  stored  energy  is  released  at  a  constant 
temperature.  The  energy  release  may  be  determined  either  as  a  function  of  time 
or  as  a  total  quantity  after  a  chosen  interval  of  annealing  time. 

(b)  Anisothermal  annealing,  where  the  temperature  is  increased  at  a  set  rate, 
and  the  released  energy  is  measured  as  a  function  of  temperature. 

Annealing  methods  have  the  added  advantage  of  providing  information  on 
the  kinetics  of  the  release  of  stored  energy. 

(2)  Reaction  method.  In  this  method  the  plastically  deformed  and  standard 
specimens  of  a  material  are  allowed  to  react  with  a  molten  substance  in  a  calori¬ 
meter.  The  difference  between  the  heats  of  reaction  of  standard  and  deformed 
specimens  is  considered  to  be  the  stored  energy  of  plastic  deformation. 

Isothermal  Annealing  Method 

Borelius  and  his  coworkers  (1953)  seem  to  have  been  the  first  to  measure 
stored  energy  by  this  method.  They  used  an  isothermal  jacket  microcalorimeter 
to  investigate  the  energy  evolved  at  100°C  by  aluminum,  copper,  and  zinc  immedi¬ 
ately  after  rolling  and  the  energy  released  at  60°C  by  aluminum  and  copper  after 
a  lapse  of  41  hours  and  45  hours  respectively.  The  calorimeter  consisted  of  a 
vapor  thermostat,  which  provided  an  isothermal  environment  for  a  chamber  con¬ 
taining  the  specimen.  The  temperature  was,  measured  by  30  junctions  of  sensitive 
differential  thermopile,  one  set  of  which  was  in  contact  with  the  chamber.  The 
other  set  was  in  contact  with  the  specimen.  In  operation,  the  sample  was  pre¬ 
heated  and  was  rapidly  placed  in  position.  The  release  of  heat  from  the  speci¬ 
men  was  indicated  by  the  temperature  difference  across  the  pile. 
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Astrom  (1956,  1964)  used  the  same  equipment  to  study  the  release  of 
stored  energy  from  aluminum  and  molybdenum  plastically  deformed  under  compression. 
Gordon  (1954,  1955)  also  used  the  Borelius  et  al .  type  of  calorimeter  and  method 
to  determine  energy  stored  and  the  kinetics  of  isothermal  release  in  the  range 
150°C  -  220°C  in  copper  specimens  deformed  under  tension.  Bailey  and  Hirsch 
(1960)  have  measured  stored  energy  in  silver  by  this  means. 

Anisothermal  Annealing  Methods 

There  are  two  basic  methods  which  may  be  used  to  obtain  results  under 
anisothermal  annealing.  Both  methods  are  normally  used  in  a  differential 
arrangement  where  the  cold  worked  sample  is  compared  to  an  annealed  sample.  In 
one  case,  the  "adiabatic  method,"  both  samples  are  heated  internally  so  as  to 
minimize  the  transfer  of  heat  between  samples  and  the  evacuated  furnace,  and 
the  stored  energy  is  given  by  the  difference  between  powers  required  to  heat 
the  two  samples.  In  the  second  method  the  heat  required  to  raise  the  tempera¬ 
ture  of  the  samples  is  determined  from  the  transfer  of  heat  from  the  furnace. 

The  two  specimens  are  placed  so  as  to  minimize  the  thermal  interaction  between 
them.  This  method  is  called  "differential  thermal  analysis." 

The  adiabatic  method  was  first  used  by  Quinney  and  Taylor  (1937)  to 
determine  stored  energies  in  copper,  70  copper  -  30  zinc  alloy,  aluminum,  iron 
and  nickel.  Welber  (1952)  and  Clarebrough  et  al .  (1952)  built  equipment  based 
on  the  design  of  Quinney  and  Taylor.  Clarebrough  and  his  coworkers  have  measured 
the  release  of  energy  stored  in  nickel  and  copper  deformed  by  various  processes. 

A  summary  of  these  investigations  has  been  given  by  Boas  (1957). 

The  differential  thermal  analysis  annealing  method  was  first  used  by 
Sato  (1931).  He  obtained  differential  heating  curves  for  cold  worked  specimens 
of  70  Cu  -  30  Zn,  silver,  copper,  aluminum,  iron,  nickel  and  brasses.  He 
annealed  standard  specimens  of  identical  materials  placed  side  by  side,  but  not 
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in  contact  with,  a  block  of  silver.  The  block  was  heated  at  the  rate  of  5°C/min 
to  the  temperature  at  which  the  standard  had  been  annealed.  The  heating  was 
repeated  without  moving  the  specimens  to  evaluate  the  corrections  necessary  to 
account  for  differences  between  the  two  specimens  not  due  to  cold  work. 

Henderson  and  Koehler  (1956)  have  used  the  technique  to  study  energy 
release  in  deformed  copper  and  brass  at  77°k.  These  materials  were  deformed 
under  compression  arid  the  energy  release  was  measured  from  77°k  to  300°k. 

Recently  Van  den  Beukel  (1963)  has  repeated  Henderson  and  Koehler's  work  and 
has  also  studied  the  energy  release  in  deformed  nickel,  silver,  and  gold. 

Reaction  Methods 

Reactions  which  have  been  used  for  measuring  stored  energy  of  plastic 
work  are  dissolution  in  aqueous  solvents  and  oxidation,  and  dissolution  in  liquid 
tin.  In  order  to  measure  the  energy  difference  between  a  deformed  and  an  unde¬ 
formed  sample  by  measuring  the  difference  between  the  heats  of  solution  in  any 
substance,  it  is  necessary  that  the  heat  of  solution  be  small,  the  solubility 
large  and  the  rate  of  solution  at  a  moderate  temperature  very  rapid.  The  only 
solvents  which  can  thus  qualify  are  liquid  metals.  The  application  of  this 
method,  though  started  in  1929  by  Smith  (1929),  has  been  applied  with  success 
starting  with  the  works  of  Bever  (1951)  and  his  associates.  The  method  they 
currently  use  is  given  by  Howlett  et  al .  (1962).  To  date,  measurements  by 
reaction  methods  have  been  limited  to  copper-gold  and  silver-gold  alloys. 

Some  other  calorimeters  which  may  be  used  to  study  release  of  deforma¬ 
tion  energy  are  those  constructed  by  Stansbury  et  al .  (1959),  Sherwin  (1964), 
Arndt  and  Fujita  (1963),  and  Calvet  and  Prat  (1963). 


2.3  Status  of  Stored  Energy  Measurements  and  Current  Problems 

The  literature  now  contains  a  considerable  number  of  investigations  in 
which  the  stored  energy  of  plastic  deformation  has  been  measured.  The  distribu¬ 
tion  of  published  papers  in  time  between  1900  and  1965  is  shown  in 
Figure  1.  However,  it  is  still  difficult  to  evaluate  critically  the  relative 
accuracy  of  the  various  methods  described  above.  There  are  at  least  three  reasons 
for  this  situation. 

(1)  Scarcity  of  data  for  any  metal  which  can  be  used  as  a  standard.  Polycrys¬ 
talline  copper  comes  closest. 

(2)  Reporting  only  of  reproducibility  and  not  accuracy  of  results  by  investi¬ 
gators.  This  is  not  adequate  since  accuracy  is  indicated  by  reproducibil ity  only 
in  the  absence  of  systematic  errors. 

(3)  The  exact  agreement  between  different  methods  should  be  possible  only  under 
certain  conditions.  These  conditions  concern  the  characteristics  of  calori¬ 
meters,  thermocouple  behavior,  thermoelastic  effects,  the  process  by  which  defor¬ 
mation  is  induced  in  a  specimen,  etc.  However,  within  limitations  imposed  by 
the  above  three  factors,  some  comparisons  are  possible;  these  have  been  made  by 
Williams  (1967)  in  a  recent  review  article. 

An  appreciation  of  the  basis  for  limited  comparisons  among  experimental 
results  may  be  obtained  by  simply  listing  variables  which  affect  the  quantity  of 
stored  energy  in  a  cold  worked  material. 

A.  Material  as  a  variable 

1.  Purity 

2.  Cdmposition  of  solid  solution  alloys 

3.  Crystallinity,  i.e.  single  crystal  or  polycrystal 

4.  If  polycrystalline,  the  effect  of  grain  size  and  grain  structure 

5.  Dislocation  characteristics,  vacancies,  etc. 

B.  Process  of  deformation  as  a  variables 

1.  Extension 

2.  Compression 
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3.  Torsion 

4.  Rolling 

5.  Wire  drawing 

6.  Cutting  processes  such  as  filing,  drilling,  etc. 

7.  Deformation  in  fatigue,  i.e.  cyclic  loading 

8.  Degree  and  rate  of  deformation 

C.  Effects  of  temperature  on  cold  work. 

Present  activity  in  the  field  is  directed  toward  evolving  a  model  or 
models  which  will  explain  the  mechanisms  of  energy  storage  and  release.  The 
present  interpretation  of  energy  storage  is  based  on  ideas  of  the  production  of 
inhomogeneous  lattice  strains  and  the  generation  of  imperfections  like  disloca¬ 
tions,  point  defects,  stacking  faults  in  hexagonal  materials,  twins  and  reduc¬ 
tion  of  order  in  alloys.  However,  a  quantitative  interpretation  of  stored 
energy  processes  remains  uncertain  primarily  due  to  lack  of  a  precise  accepted 
definition  of  cold  work  and  to  the  degree  of  inaccuracy  in  measurement  of 
stored  energy. 

It  will  be  seen  from  the  above  .description  that  the  stored  energy  of 
cold  work  has  been  measured  exclusively  in  metals.  Even  among  metals,  copper, 
nickel,  aluminum,  silver  and  iron,  have  been  most  commonly  investigated; 
infrequent  or  solitary  measurements  of  stored  energy  of  cold  work  in  zinc,  lead 
tin  do  exist.  The  stored  energy  of  cold  work  has  also  been  investigated  in  a 
few  alloys,  notably  in  Cu-Ag,  Cu-Au,  and  Cu-Ni.  The  absence  of  any  such  inves¬ 
tigation  in  ionic  compounds  is  striking  considering  the  large  amount  of  work 
done  in  LiF,  MgO  and  other  ionic  solids  (see  for  example,  (Gilman,  1969)). 

Laboratories  where  investigations  on  the  stored  energy  of  cold  work  and 
its  release  mechanism  are  being  carried  out  at  present  are  few.  The  prominent 
laboratories  and  persons  involved  in  these  measurements  are: 

(1)  Oak  Ridge  National  Laboratory,  Oak  Ridge,  Tenn.  R.  0.  Williams,  A. 

Wolfenden,  and  R.  A.  Vandermeer. 
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(2)  Dept,  of  Metallurgy,  MIT,  Cambridge,  Mass.  M.  B.  Bever  and  A.  L.  Titchener 

(3)  Dept,  of  Metallurgical  Engineering,  I IT ,  Illinois,  P.  Gordon 

(4)  Cavendish  Laboratory,  Cambridge,  U.K. ,  J.  E.  Bailey  and  P.  B.  Hirsch 

(5)  Dept,  of  Metallurgy,  Univ.  of  Liverpool,  U.K.,  A.  S.  Appleton 

(6)  Div.  of  Tribophysics,  C.S.I.R.O.,  Univ.  of  Melbourne,  Australia, 

L.  M.  Clarebrough,  M.  E.  Hargreaves,  M.  H.  Loretto  and  A.  K.  Head 

Present  ideas  about  the  mechanisms  of  storage  of  cold  work  energy  and 
its  release  are  given  in  what  follows. 

The  small  percentage  (5  to  10%)  of  the  energy  expended  in  plastically 
deforming  a  metal,  which  is  stored  in  the  deformed  metal,  is  associated  with 
defects  generated  during  the  deformation.  A  significant  part  of  the  stored 
energy  is  associated  with  the  distribution  and  density  of  dislocations  in  the 
bulk  of  the  deformed  metal.  The  remaining  part  of  the  stored  energy  is  associ¬ 
ated  with  imperfections  like  point  defects,  vacancies,  twins,  stacking  faults, 
etc.  However,  this  partition  of  stored  energy  among  the  various  imperfections 
as  well  as  the  absolute  value  of  the  stored  energy  is  dependent  on  purity, 
degree  of  deformation,  temperature  of  deformation,  and  grain  size.  In  general, 
the  energy  stored  during  plastic  deformation  of  a  metal  increases  with  the  degree 
of  deformation,  and  may  or  may  not  reach  a  saturation  before  fracturing.  Also 
for  a  given  extent  of  deformation,  the  value  of  stored  energy  increases  with 
decreasing  temperature  of  deformation  of  a  metal.  In  general,  the  higher  the 
purity  of  a  metal,  the  lower  the  stored  energy  of. cold  work.  And  lastly,  for 
metals,  stored  energy  for  a  given  deformation  increases  with  decrease  in  grain 
size  for  small  strains  (<0.20)  and  then  becomes  relatively  independent  of  grain 
size  for  large  strains. 

The  inference  about  apportionment  of  stored  energy  among  the  various 
imperfections  in  a  metal  is  made  from  the  correlation  of  release  of  stored 
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energy  with  changes  in  other  physical  properties,  e.g.,  changes  in  electrical 
resistivity,  density  and  with  structural  studies  using  conventional  metallo- 
graphic  techniques  and  electron  transmission  microscopy. 

The  release  of  stored  energy  proceeds  in  two  or  more  stages.  However  the 
various  processes  by  which  the  annealed  state  is  reached  are  classified  into  two 
categories  --  recovery  and  recrystallization.  The  apportionment  of  the  released 
energy  into  these  two  classes  is  based  on  the  definition  that  recovery  encompasses 
all  those  annealing  phenomena  which  precede  recrystallization  in  any  local  region 
of  the  material.  On  the  other  hand,  recrystallization  is  defined  as  the  appear¬ 
ance  of  comparatively  strain-free  grains  which  are  separated  from  the  deformed 
matrix  or  from  other  recrystallized  grains  by  high-angle  grain  boundaries  and 
alternatively,  recrystallization  may  be  defined  as  the  appearance  of  virtually 
strain-free  grains  which  have  grown  in  the  deformed  matrix  by  the  movement  of 
large  angle  grain  boundaries. 

A  common  feature  in  all  the  investigations  on  metals  and  alloys  is  that 
recrystallization  is  accompanied  by  a  marked  release  of  stored  energy  as  well  as 
-a  decrease  in  electrical  resistivity  and  an  increase  in  density.  Electron  micro¬ 
scope  observations  indicate  that  this  release  of  energy  and  concomittant  changes 
in  the  electrical  resistivity  and  density  are  associated  with  a  change  in  the 
distribution  and  density  of  dislocations.  Hence  the  magnitude  of  energy  released 
during  recrystallization  can  be  associated  with  the  energy  of  the  dislocations. 
Bailey  and  Hirsch  (1960)  seem  to  have  been  the  first  to  establish  such  an 
association,  which  they  did  for  polycrystalline  silver.  However,  Clarebrough 
et  al .  (1961)  disputed  the  conclusion  of  Bailey  and  Hirsch  that  theory  of  flow 
stress  involving  piled  up  groups  of  dislocations  is  not  applicable  to  silver. 

For  recovery  the  experimental  situation  is  not  so  well  defined.  The 
ratio  of  stored  energy  released  during  recovery  to  total  stored  energy  varies 
between  0.03  to  0.7  depending  on  the  purity  of  material,  strain,  and  temperature 
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of  deformation.  Evaluations  of  changes  in  other  properties  accompanying  recovery 
in  metals  deformed  at  room  temperature  have  been  aided  insignificantly  by 
structural  observations.  For  metals  deformed  at  the  temperature  of  liquid  nitro¬ 
gen  the  energy  stored  is  four  times  that  for  deformation  at  room  temperature. 

The  annealing  of  point  defects  is  assumed  to  contribute  to  the  large  releases 
of  energy  associated  with  recovery  below  room  temperature,  in  addition  to  the 
contribution  of  other  imperfections.  Bailey  and  Hirsch  (1960)  could  not  detect 
any  change  in  the  dislocation  density  and  distribution  before  and  after  recovery 
in  silver,  although  during  the  recovery  half  the  total  stored  energy  was  liberated. 

In  alloys  a  major  source  of  energy  release  during  recovery  may  be 
associated  with  the  return  of  short  range  order  destroyed  by  deformation. 

To  sum  up,  existing  experimental  measurements  on  the  stored  energy  of 
cold  work  and  the  release  of  the  energy  with  concomittant  variation  in  the  physi¬ 
cal  structural  properties  in  the  most  extensively  studied  metal,  i.e.  copper,  is 
not  sufficient  to  test  validity  of  the  two  existing  theories  of  work  hardening 
of  pure  metals.  In  one  theory,  due  to  Friedel  (1955)  and  Seeger  et  al .  (1957), 
the  flow  stress  is  determined  by  long-range  stresses  due  to  groups  of  piled  up 
dislocations.  In  the  second  theory,  due  to  Cottrell  (1953),  Basinski  (1959), 
and  Hirsch  (1958),  the  flow  stress  is  determined  by  short-range  interactions 
between  crossing  dislocations. 
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III.  PLASTIC  DEFORMATION  IN  UNIAXIAL  STRAIN 

3.1  Definitions  and  Basic  Relations 

The  uniaxial  strain  condition  applies  to  plane  wave  compression  and 
rarefaction.  If  wave  propagation  is  in  the  x-direction,  the  work  of  uniaxial 
compression  is 

dW  =  -Vpx  dex  =  -px  dV  (1) 

where  p  =  -o  =  compressive  stress  in  the  x-direction  and  de  =  dV/V.  Prin- 
cipal  axes  of  stress  and  strain  lie  in  and  normal  to  the  direction  of  propagation. 
If  E  is  internal  energy  of  unit  mass, 

dE  =  dQ  -  Vpx  dex  (2) 


where  dQ  is  heat  communicated  to  the  mass  element  by  conduction  or  other  means. 
The  work  term  of  Eq.  (1)  can  be  expressed  as  work  of  compression  plus  work  of 
deformation  by  introducing  stress  and  strain  deviators,  S-  and  e • : 

J  J 


dW  = 

-Vp  de  +  V  l  S.  de. 

j  0  J 

(3) 

where 

sj  ■ 

-Pj  +  P  ) 

(4a) 

)  j  *  x,  y,  z. 

Ej  ‘  9/3  * 

(4b) 

s*  ■ 

4t/3,  S  »  Sz  =  -2t/3 

(5a) 

T  - 

-(Px  -  Py)/2 

(5b) 

ex  = 

V3>  ey  '  ez  -  -Ex/3-  de  '  dex- 

(5c) 

Using  these  relations,  Eq.  (3)  becomes 
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dW  =  -Vpdex  +  V[(4t/3)  •  (2/3)dex  -  2(2t/3)  (-de/3)] 


=  -Vpdex 

+  4Vxde  /3 

A 

(6) 

=  -Vpdex 

(1  -4x/3p) 

(7) 

=  -Vpdex 

0  -  R„) 

(8) 

where  RD  is  the  ratio  of  work  of  deformation  to  compressive  work.  For  elastic 
changes  the  constitutive  relations  can  be  written  in  incremental  form  to  account 
for  changes  in  the  elastic  moduli  with  strain: 


dp 

=  -Kde  =  -Kde 

A 

(9) 

dsj 

=  2vde3 

(10) 

dox 

=  (x  +  2y)dex 

(11) 

d°y 

=  XdEx  »  dc2 

(12) 

dx 

=  pdex 

(13) 

where  K  is  bulk  modulus  and  x,  y  are  the  Lame  constants.  Values  of  K,  A,  y 

depend  on  whether  strain  occurs  at  constant  temperature  or  constant  entropy. 

ITT  S  S  S 

Isothermal  moduli  will  be  denoted  K  ,  X  ,  y  ;  isentropic  moduli  K  ,  x  ,  y  . 

When  the  superscript  is  omitted  the  meaning  is  obvious  or  the  question  is 

indifferent. 

We  can  invoke  the  condition  of  uniaxial  strain  to  obtain: 


Q- 

TD 

<< 

il 

-doy  =  vdpx/(l-v)  =  dp2 

04) 

dp  = 

(dpx/3)  (l+v)/(l-v) 

(15) 

dx  = 

-(dpx/2)  (l-2v)/(l-v) 

(16) 

where  v  =  (x/2)/(x+y)  =  Poisson's  ratio,  which  may  depend  on  strain  and  tem¬ 


perature  or  entropy. 
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For  small  strains,  Eqs.  (11)  -  (13)  can  be  integrated  assuming  moduli  to 


be  constant: 

ax  =  (X  +  2u)  ex 


(17) 


(18) 


T  =  yex 

ex  =  (V-V0)/V0 


(19) 

(20) 


where  V  is  specific  volume  and  V'  denotes  a  reference,  stress  and  strain-free 
state.  For  this  case 

Rd  =  2(1  -  2v)/(l  +v)  (21) 

(see  Eq.  (8)).  The  dependence  of  Rp  on  v  is  shown  in  Fig.  2. 

This  discussion  of  elastic  strain  provides  a  framework  within  which 
questions  of  plasticity  may  be  discussed  and  it  is  important  to  note  that  the 
work  of  plastic  deformation  is  bounded  by  the  work  of  volume  compression  and 
that  of  elastic  deformation.  This  is  illustrated  in  Fig.  3  where  curves  of 
hydrostatic  and  uniaxial  compression  under  adiabatic  conditions  are  compared. 

For  elastic  strain,  adiabatic  conditions  are  also  isentropic,  so  the  slopes  of 
0M  and  ON  are  -(V  /V)K^  and  (-V  /V)  (x^  +  2y^),  respectively.  Since  p  = 

U  U  A 

p  -  4t/3,  the  work  of  deformation  is 


Wd  =  -/  PxdV  +  JpdV  =  (4/3)  / t  dV  (22) 

which  is  just  VQ  times  the  area  of  the  sliver  between  ON  and  0M  of  Fig.  3.  As 
noted  in  Fig.  2,  this  may  be  small  or  large  depending  on  the  value  of  v. 

When  the  strain  cycle  is  completed  by  allowing  the  material  to  expand 
back  to  p  or  p  =  0,  the  elastically  strained  sample  under  adiabatic  conditions 

A 

simply  returns  to  the  point  0  along  the  curve  of  compression,  OA  or  OB.  The 
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Figure  3  Adiabatic  elastic  compression  of  an  isotropic  solid 
in  uniaxial  strain  and  hydrostatical ly . 
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plastically  deforming  material,  on  the  other  hand,  completes  part  of  a  hysteresis 
loop  as  shown  in  Fig.  4.  The  sections  of  this  loop  are  formed  according  to  the 
following  process.  Starting  at  point  0  the  sample  is  compressed  uniaxial ly  until 
it  starts  to  deform  plastically  at  A.  Along  AB  the  plastic  deformation  is  con¬ 
tinuous  until  the  maximum  value  of  px  is  reached.  At  B  the  applied  stress  is 
allowed  to  start  decreasing,  a  -  a  is  reduced  and  expansion  occurs  elastically. 

x  y 

At  C,  a  -  a  has  once  more  reached  the  yield  value  and  plastic  deformation  con- 
y  x 

tinues  from  C  to  D.  There  are  now  two  questions  to  be  considered: 

(i)  What  is  the  amount  of  plastic  work  done  along  the  sections  AB 

and  CD?  (Ql) 

(ii)  Can  a  thermodynamic  state  be  defined  at  each  point?  (Q2) 

We  consider  these  questions  in  the  order  stated. 
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HYDROSTATIC 


Figure  4  Adiabatic  compression  and  expansion  of  an  elastic- 
plastic  solid  in  uniaxial  strain,  OABCD. 
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3.2  Computation  of  the  Plastic  Work 

An  important  convention  in  the  theory  of  plasticity  is  the  assumption 
that  when  plastic  and  elastic  strain  are  occurring  simultaneously,  each  increment 
in  total  strain  can  be  decomposed  into  a  sum  of  plastic  and  elastic  increments, 
i.e. 

de  =  dee  +  de^  (23) 

where  superscripts  "e"  and  "p"  stand  for  "elastic"  and  "plastic." 

This  convention  is  discussed  in  Appendix  A.  It  is  also  common  to  assume 
that  there  is  no  plastic  contribution  to  the  density  change,  so  that 


de  .  ' =  de  .e  +  de  =  de .  +  de/3 

J  J  J  J 


de  =  de'  ; 


de •  =  de •  +  de . 

J  J  J 


=  0 


(24a) 

(24b) 


Then  increments  in  the  work  of  plastic  and  elastic  deformation  are  given, 
respectively,  by 


dwdp  ■  U  so  de/ 

plastic 

(25) 

dWde  *  nSlde,e 

J  J  J 

elastic 

(26) 

Since  dW^  is  the  only  plastic  work  in  this  process,  i.e.  volume  compression  is 
purely  elastic,  we  write  henceforth 

dW_  =  dW  ,  . 

P  dp 

We  also  assume  that  stresses  are  supported  only  by  the  elastic  strains, 


dS,  =  2y  de.. 

J  J 


(27) 


A  physical  basis  for  Eq.  (27)  is  discussed  in  Appendix  B. 
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The  condition  of  uniaxial  strain  requires  only  that  total  strain  ey 
and  e2  vanish,  i.e. 

deyP  +  deye  =  0  (28a) 

dezP  +  deze  =  0  .  (28b) 

Because  of  symmetry  about  the  x-axis 

d£ye  =  deze  .  (29) 

From  Eqs.  (5a),  (10)  and  (27) 


dSx 

=  4dx/3 

2y  dexe 

dsy 

=  -2dt/3  = 

2y  deye 

dSx 

-  dSy  =  2dx 

-  2y  d(exe 

dr 

=  2ydye 

where  dye  =  d(exe  -  eye)/2  =  (dex6  ’  dey6^2' 

From  Eq.  (6)  the  total  incremental  work  of  deformation  is 


(30) 

(31) 


(32) 

(33) 


dWj  =  4Vt  dex/3.  (34) 

From  Eqs.  (26)  and  (30)  -  (32)  the  incremental  work  of  elastic  deformation  is 
dWde  =  8Vt  dye/3 

=  (4Vx/3y)dT  .  (35) 


The  incremental  work  of  plastic  deformation  is  obtained  by  subtracting  Eq. 

(35)  from  (34): 

dW  =  (4/3)  Vt  (de  -  d-r/y)  (36) 

P  x 
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where  t  =  ±Y/2  depending  on  whether  the  process  is  compressive  or  expansive. 

Y  >  0  is  the  yield  stress  as  given  by  the  von  Mises  condition. 

Eq.  (36)  can  be  obtained  directly  from  Eqs.  (4b),  (24a)  and  (25).  It  can  also  be 
written  in  the  form  of  Eq.  (C5)  in  Appendix. C: 

dW  =  2Vx  de  p  '  (37) 

P  X 

If  y*5  =  (e  p  -  e  p)/2,  Eq.  (37)  can  also  be  written  as 
x  y 

dWp  =  (8/3)Vx  d  yP  (38) 

The  increment  in  elastic  work  represented  by  Eq.  (35)  can  easily  be  seen  to  be 
reversible: 

dWde  =  (2V/3y)  d(x2).  (39) 

As  | t )  goes  up  and  down,  so  does  Wde. 

The  plastic  work  is  irreversible  by  definition, 

dWp  >0.  (40) 

This  is  readily  illustrated  with  reference  to  Fig.  4.  During  elastic  compression 
along  OA  in  Fig.  4,  de  p  =  0,  t  <  0,  dr  <  0,  dW  =  0.  Along  AB,  t  <0  and 

X  p 

de  <  0.  For  a  material  which  doesn't  work-harden,  dx  =  0,  so  dW  >  0.  If 
x  P 

work-hardening  occurs,  dx  <  0,  but  dx/y  cannot  exceed  de  in  magnitude  since 

<  A 

de  =  dx/y  when  compression  is  elastic.  So  in  that  case,  too,  dW  >  0.  Along 
x  p 

BC  the  material  is  expanding  elastically,  x  is  changing  from  negative  to  positive 
sign,  dexP  =  0  and  dWp  =  0.  Along  CD,  x  >  0,  dex  >  0  and  again  dWp  >  0  from 
Eq.  (36)  by  the  same  argument  applied  to  the  path  AB.  Also  de  p  >  0  along  CD. 

A 

Eq.  (36)  or  (37)  provides  an  answer  to  question  T,  raised  above.  We 
now  turn  attention  to  Q2. 
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IV.  THERMODYNAMIC  STATES  IN  PLASTIC  FLOW 


The  simplest  observations  of  plastic  yield  tell  us  a  great  deal  about  the 
process.  Bend  a  piece  of  iron  wire  back  and  forth  between  your  fingers  a  few 
times.  It  gets  hot  and  is  still  solid  iron.  The  first  result  is  obvious;  the 
second  is  verified  by  chemical  and  metallurgical  analysis.  It  is,  in  fact,  very 
hard  to  detect  any  changes  in  the  material,  but  we  know  that  some  have  occurred 
since  repeated  bending  breaks  the  wire,  yet  it  does  not  break  at  the  first  bending. 
Property  changes  due  to  bending  are  more  obvious  in  a  piece  of  copper.  Bending 
produces  marked  increases  in  bending  resistance  before  it  breaks. 

In  developing  a  thermodynamic  theory  to  deal  with  such  conditions  we  can 
assume  either  that  plastic  work  changes  the  material  or  that  it  does  not.  It  is 
clear  from  our  simple  experiment  that  the  latter  assumption  is  a  fairly  good  one 
and  we  accept  it  for  the  moment  and  consider  its  consequences. 

4.1  Plastic  Work  Converted  Entirely  to  Heat 

Even  though  this  be  true,  there  is  a  difficulty  associated  with  plastic 
strain  which  is  illustrated  in  Fig.  4.  On  reversal  of  stress  at  any  point,  say 
B,  the  stress-strain  curve  is  not  retraced,  even  though  the  process  be  isothermal. 
Yet  if  we  persist  in  the  reversal  we  can  find  a  sequence  of  applied  forces  which 
will  return  the  material  to  its  initial  state  under  this  condition  that  plastic 
work  produces  no  material  changes.  This  is  illustrated  in  Figs.  6  and  7,  to  be 
discussed  farther  on  in  this  section.  This  is  a  kind  of  geometric  irreversibility 
to  be  distinguished  from  irreversibility  of  material  properties.  It  complicates 
description  of  thermodynamic  processes  but  does  not  essentially  change  the  postulate 
that  plastic  working  has  put  the  material  in  a  thermodynamic  state  characterized 
by  the  added  heat,  aQ.  The  situation  is  somewhat  like  that  of  a  fluid.  If  we 

deform  a  viscous  fluid  at  constant  density,  we  recognize  that  heating  occurs, 
but  we  suppose  the  final  state  to  be  describable  in  terms  of  thermodynamic 
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parameters  of  the  original  fluid.  The  fact  that  particular  molecules  have 

1 

different  neighbors  is  considered  inconsequential. 

In  the  same  way  it  is  inconsequential  that  neighbors  may  have  changed 
after  plastic  deformation;  thermodynamic  parameters  are  the  same.  It  is,  in  a 
sense,  unimportant  that  the  macroscopic  form  of  the  material  has  changed  and 
that  this  makes  a  difference  in  the  external  behavior.  This  difference  is 
observable  in  the  solid  and  not  in  the  fluid  because. the  solid  atoms  undergo 
elastic  strain  during  deformation,  whereas  the  fluid  atoms  do  not. 

With  such  considerations  as  foundation,  we  are  led  to  a  Fundamental 
Assumption: 

There  is  always  an  elastically  strained  state 
imbedded  in  every  deformed  state.  This  elastic 
state  is  inherently  reversible  and  can  be  recovered 
by  suitable  processes. 

This  assumption  implies  that  an  internal  energy,  E,  exists  and  that  it 
is  a  function  of  the  elastic  strains,  ee,  entropy,  S,  and  perhaps  of  inter 
nal  or  microscopic  strain  variables.  A  representation  of  the  internal  energy 
function  can  be  obtained  by  noting  that  each  increment  of  work  on  the  system 
is  the  sum  of  an  elastic  and  a  plastic  increment.  By  the  Fundamental  Assump¬ 
tion  the  elastic  increment  is  recoverable,  so  increments  in  internal  energy 
can  be  written  as 

dE  =  TdS  +  dWe  (41) 

where  dWg  depends  on  the  elastic  strain  increments.  According  to  the  first 
law  this  increase  in  internal  energy  must  equal  the  total  work  done  on  the 
system  plus  the  heat  absorbed  from  the  surroundings: 
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dW  +  dQ  . 


(42) 


TdS  +  dWQ  = 
e 

This  leads  directly  to  the  Gibbs  relation  for  the  increase  in  entropy  pro¬ 
duced  in  the  process: 

dS  =  (dWp  +  dQ)  /  T  (43) 

where  dW„  =  dW  -  did  . 

P  e 

According  to  Eq.  (41)  and  the  foregoing  assumptions,  internal  energy  can 
be  expressed  in  terms  of  entropy  and  the  elastic  strains.  For  example 

E  ■  E(S,  Eye) 

Changes  in  entropy  result  solely  from  the  agitations  and  elastic  displacement 
of  atoms;  there  is  no  entropy  directly  associated  with  the  plastic  strain. 
Nonetheless  an  apparent  dependence  of  state  on  plastic  strains  may  exist 
because  of  the  procedure  required  to  separate  dW  into  elastic  and  plastic 
parts.  For  example,  from  Eqs.  (1)  and  (37)  we  find  that 

dW  =  dW  -  dWn  =  Vo  de  -  2Vxde  P 

c  p  A  X  X 

so  that 

E  =  E(S,  ex,  e/). 

The  appearance  of  exP  in  this  potential  is  illusory.  It  does  not  imply  a 
physical  dependence  of  E  on  plastic  strain.  Failure  to  realize  this  fact  has 
led  to  some  confusion  in  other  work. 

It  will  be  useful  to  tabulate  various  forms  of  the  internal  energy,  but 

\ 

first  we  note  that  increments  in  the  Helmholtz  potential  will  be  in  the  form 

dA  =  -SdT  +  dWe  (44) 
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i 

and  that  the  entropy  and  temperature  are  given  by  ; 

S'  =  -(3A/3T)  p 

(45) 

£ 

T  =  (3E/3S)  p 

(46) 

e 

where  subscript  ee  indicates  that  all  elastic  strains  are  to  be  held 

constant. 

Using  the  equations  developed  earlier,  we  have  the  following  representations: 

From  Eqs.  (6)  and  (35) , 

dW  =  -Vpde  +  8VxdYe/3 

ti  A 

(47) 

E  =  E(S ,  ex,  ye) 

(48) 

A  =  A(T ,  ex>  ye) . 

(49) 

Since  VdE  =  dV,  this  is  equivalent  to  a  representation  in  which 

A 

V  and 

ye  are  the  mechanical  variables. 

From  Eqs.  (1 )  and  (37) , 

dW  =  dW  -  dW  =  Va  deY  -  2VxdEvP 

c  [J  A  A  A 

(50) 

E  ■  E(S,  ex,  c/) 

(51) 

A  -  A(T.  ex,  c/). 

(52) 

From  Eqs.  (1 )  and  (38) 

dW  =  dW  -  dW  =  Va  de  -  8VxdYP/3 
s  p  XX 

(53) 

> 

E  =  E(S,  ex,  yP) 

(54) 

A  =  A(T ,  ex,  yP)- 

(55) 
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A  representation  which  shows  the  dependence  of  state  functions  on  elastic  strains 
and  temperature  exclusively  is 

dWe  =  V(ax  dexe  +  2 ay  dEye)  (56) 

so  that 

E  -  E(S,  exe,  Eye)  (57) 

A  =  A(T ,  Exe,  sye)  (58) 

Since  entropy  is  defined  as 

S  =  - 3A/ 3T , 

it  depends  on  the  same  variables  as  A. 

The  functions  E,  A,  S,  etc.  are  state  functions,  so  their  differentials  are 
exact  and  lead  to  various  Maxwell -like  relations  and  compatibility  conditions 
which  are  later  shown  to  be  useful. 

Some  common  thermodynamic  observables  take  unusual  forms  when  these 
representations  are  used.  Specific  heat  is  an  example. 

From  Eq.  (49) , 

C„  =  T(3S/3T)  =  T( 3S/3T)  0  +  T(3S/3ye)T  (3Ye/3T)c  (59a) 

X  EX,Y  x  r  X 

From  Eq.  (52) , 

Cw  =  T( 3S/3T)  =  T(3S/3T)  +  T(3S/3e  P)  T  ( 3e  P/3T)  (59b) 

V  £x  ex,ExP  x  ex’’  x  X 

From  Eq.  (55) , 

Cu  =  T(3S/3T)  =  T(3S/3T)  _  +  T(3S/3yP)  t  (3yP/3T)  (59c) 

v  £X  ex,YP  £X,  £X 

For  Eq.  (56)  C y  is  not  defined  but  can  be  obtained  by  invoking  the  condition 

de  =  de  e  +  2de  e,  (60) 

x  y 
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which  follows  directly  from  Eq.  (24a)  and  the  symmetry  of  uniaxial  strain.  In 

this  way  one  of  the  three  forms  of  Eq.  (59)  will  be  obtained. 

One  may  reasonably  ask  whether  derivatives  like  (3S/3e  P)T  which 

x  I  *exs 

appears  in  Eq.  (59b),  are  physically  meaningful  in  view  of  the  nature 
of  plastic  strain.  Since  plastic  strain  enters  in  a  contrived  way,  as  explained 
in  the  paragraph  following  Eq.  (43),  they  can  always  be  expressed  in  terms  of 
more  obviously  thermodynamic  derivatives.  For  example,  calculate  dS  from 
Eqs.  (52)  and  (58): 

dS  =  (3S/3T)  dT  +  (3S/9ex)  dex 

ex’ex  T’ex 

+  (9S/3exP)T  dexP 

X 

dS  =  (3S/3T)  e  e  dT  +  (3S/3exe)  c le / 

Ex  ,ey  T,Ey 

+  (3S/3T)  de  e  . 

TjEx 

Since  we  wish  to  obtain  3S/3e  p  in  terms  of  elastic  strain  derivatives,  write 

A 

de  e  and  de..e  in  terms  of  de  and  de  p;  i.e. 
x  y  xx 

de  ^  =  de  -  de  P 

X  X  X 

dey6  =  -deyP  =  dexP/2  . 

•  Then  equating  the  two  values  of  dS  gives 
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0  =  [(3S/3T)  -  (3S/8T)  e  e]dT  +  [(3S/3ex) 

Vx  Ex  ey  T'£x 

-(3S/3Exe)  e]dex  *  [<*/..„% 

T’Ey  x 

-(l/2)(3S/3c  e)  +  (»S/3Exe)  e]  dsxp  .  (61) 

t-ex  T,ey 

Specification  of  the  deformation  process  requires  specification  of  dT  and  two 
strain  variables,  so  coefficients  of  dT,  dexe,  dexP  in  Eq.  (61)  must  vanish 
independently;  therefore 


(3S/3cxP)T  e  =  (1/2)  (as/e /)  -  (3S/3e  e) 

x  1  ,ex  y  T  »ew  T,< 


(62) 


From  Eqs.  (44)  and  (56) 


dA  =  -SdT  +  Voxdexe  +  2Vaydeye 


S  =  -9A/9T;  Vax  =  9A/9ex 


9S/9e: 


-  92A/3T9e  e  =  -  (9Vax/3T)  e 

ex  ey 


9S/ 3e 


2(3Vo  /9T)  e  e 
ey  £x 


Then 


(9S/9exP)T>( 


V(9ax/9T)  e  e  -  V(9ay/9T)  e 
ex  ey  ex  Gy 


=  2V(9t/9T) 


e  e 


(63) 


The  derivative  on  the  r.h.s.  of  Eq.  (63)  is  taken  with  e  e,e  e  constant.  Since 

x  y 

this  insures  that  no  strain  is  occurring, 

( 3t/3T)  e  e  =  (3t/3T)  e  =  (3t/3T)  e  etc. 

Ex  ey  Vex  VY 

In  order  to  evaluate  changes  in  entropy  under  the  assumptions  of  this 
section,  dQ  of  Eq.  (43)  must  be  specified.  For  an  adiabatic  process,  dQ  =  0, 
but  the  sample  gets  hotter  as  the  process  continues,  so  a  repeating  cycle 
doesn't  exist. 

To  calculate  dQ  for  an  isothermal  process  consider  incremental  com¬ 
pression  from  a  reference  state  under  three  different  conditions.  Fig.  5. 
Compression  along  AD  is  adiabatic  and  is  accompanied  by  plastic  work  aW  ; 

r 

no  heat  is  absorbed.  Compression  along  AC  is  isentropic;  the  amount  of  heat 
removed  just  compensates  for  the  plastic  work.  Compression  along  AB  is  iso¬ 
thermal  and  it  is  along  this  path  that  we  wish  to  calculate  aQ.  For  the 
situation  shown  aQ  <  0,  i.e.  heat  must  be  removed  to  compensate  for  the 
plastic  work.  By  the  first  law  the  change  in  internal  energy  in  going  from  A 
to  B  is  equal  to  the  total  work  done  plus  the  heat  absorbed, 

EB  -  EA  =  dW  +  dQ  (64) 

where 

AW  =  Vox  dex 

According  to  Appendix  F,  dQ  can  be  expressed  in  terms  of  incremental  changes 
in  specific  volume,  elastic  strain  and  plastic  work: 

dQ  =  T(3p/9T)  dV  -  (8VT/3)(9t/9T)  e  dYe  -  dWn  (65) 

V, ye  V,Ye  P 
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V+dV 


V 


Figure  5  Heat  absorbed  in  isothermal  compression 


1 


The  net  entropy  production  in  the  sample  is  obtained  by  combining  Eqs.  (43) 
and  (65): 

dS  =  (l/T)(dW  +  dQ) 

=  (3p/3T)  p  dV  -  (8V/3)(3x/3T)  ,  p  dTe  (66) 

This  is  a  remarkable  result  and  one  which  has  been  noted  by  Bridgman  for 
uniaxial  stress:  Entropy  absorbed  by  the  sample  under  isothermal  conditions 
is  given  by  the  elastic  changes  alone,  even  in  the  presence  of  plastic  work. 

Equation  (66)  may  be  obtained  more  directly  with  minimal  reference  to 
the  physical  processes  by  using  the  fact  that  entropy  is  a  state  function,  so 
changes  in  entropy  can  be  obtained  by  direct  differentiation.  Thus  from  Eq. 
(49)  with  Vde  replaced  by  dV, 

A 

S  =  S(T,  V,  Ye) 

dS  =  '(3S/3T)  e  dT  +  (3S/3V)  dV  +  (3S/3Ye)T  v  dY6  (67) 

V,Y  T»y 

Since  S  is  derived  from  the  Helmholtz  potential  A  and 

dA  =  -SdT  -  pdV  +  (8Vt/3)  dYe, 

the  last  two  coefficients  in  Eq.  (67)  can  be  obtained  by  equating  cross¬ 
derivatives  of  A: 

3S/3V  =  -32A/ ST3V  =  3p/3T  (68) 

aS/3Ye  =  -a2A/3T3Ye  =  - ( 8V/3 ) ( 3t/  )  (69) 

Setting  dT  =  0  and  substituting  Eqs.  (68)  and  (69)  into  (67)  yields  Eq.  (66) 
While  this  procedure  is  economical  in  mathematical  operations,  it  does  not 
reveal  the  physical  content  of  the  result  so  clearly  as  the  earlier  procedure. 
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Since  dS  is  an  exact  differential,  the  cross  derivatives  of  its  coeffi¬ 
cients  must  be  equal.  Thus 

(acv  /3V)  e  =  (T  32p/3T2)  e  '  (70a) 

T,y  V,y 

(9CVy/3ye)VJ  =  -(8V/3)(a2T/8T2)v>ye  (70b) 

32p/8T9ye  =  -(8/3)  3(V3t/3T)/3 V  (70c) 

The  bulk  modulus,  K,  may  depend  on  ye,  thus  leading  to  a  dependence  of  p  on 
ye,  but  this  is  likely  to  be  small.  Setting  it  equal  to  zero  implies  that 

3(V3t/3T)/3 V  =  0 

or 

V(3x/3T)  =  f(ye,T)  (71) 

V,y6 

In  Eqs.  (70)  C„  has  been  written  for  (1/T)(3S/3T)  Q.  Here  again  the 

y  e  V,Y 
dependence  of  Cyy  on  y  is  likely  to  be  small  so 

(32t/3T2)  «  0.  (72) 

V,ye 

Most  materials  will  be  found  to  yield  before  the  linear  relation  between  t 
and  y6  is  violated.  Setting  y  =  2y  ye  with  y  independent  of  y6,  Eq.  (71) 
becomes 

(3y/3T)  p  =  g(T)/V  (73) 

V,ye 

while  Eq.  (72)  suggests  that 

(32u/3T2)  0  -  0.  (74) 

V,ye 
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Equations  (73)  and  (74)  together  imply  that 

9y/9T  =  b/V  (75a) 

or 

y  =  y0  +  b(T/V  -  Tq/Vq)  +  4>(V)  (75b) 

where  ( Vq)  =  0.  If  experiments  reveal  that  9y/9T  depends  on  temperature 

(which  it  certainly  does  near  the  melting  temperature)  then  the  assumption  that 
C is  independent  of  ye  is  invalid. 

For  computational  purposes  it  is  convenient  to  write 

(9P/9T)  e  =  rCv  /V  (76) 

VjY  •  ' 

where  r  is  the  Gruneisen  parameter.  This  coefficient  is  slowly  varying  and 
normally,  but  not  necessarily,  positive.  Taking  x  3  2yye  as  before,  the 
second  coefficient  in  Eq.  (66)  becomes 

-  (16V  y6/3)  (9y/9T)  (77) 

V,y6 

Since  9y/9T  is  normally  negative,  this  coefficient  has  the  sign  of  ye ,  which 
is  frequently,  but  not  always,  negative  in  compression  and  positive  in  tension. 
Taking  rCy  and  V  9y/9T  to  be  constants,  we  have  for  the  entropy  change 

<f>  dS  =  TCV  $dex  +  (8V/3)  1 9y/9T |  $d(Y®)2  =  0  (78) 

This  illustrates  in  a  particular  case  the  general  result  which  follows  from 
the  assumption  that  E  and  A  are  state  functions. 

Equation  (78)  refers  to  entropy  change  in  the  sample  during  an  iso¬ 
thermal  cycle.  The  entropy  of  the  reservoir  in  which  the  sample  is  imbedded 
will  be  found  to  increase  in  a  cycle  in  the  amount 

AS  =  (dWp/T) . 
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A  simple  example  of  the  calculation  outlined  in  this  section  is  detailed 
in  Appendix  D  and  in  Figs.  6  and  7.  In  this  example  it  is  assumed  that  the 
material  is  linearly  elastic  up  to  the  yield  point  (point  A  in  Fig.  4)  and 
that  the  yield  stress,  Y  >  0,  is  constant  throughout  the  isothermal  process. 
Under  these  conditions  the  cycle  0ABF0  of  Fig.  4  can  be  followed  explicitly 
but  simply.  The  turning  point  F  is  determined  by  equating  the  tensile  stress 
developed  along  the  curve  of  elastic  expansion,  OF,  to  the  stress  along  DF. 

With  F  determined  in  this  way,  stress,  elastic  strain,  plastic  strain  and 
entropy  of  the  sample  return  to  their  starting  values  at  0,  whereas  the  plastic 
work  of  deformation  and  the  entropy  of  the  heat  reservoir  do  not.  These 
results  are  shown  in  Figs.  6  and  7.  With  non-linear  elastic  behavior  the 
computation  is  more  difficult  but  results  are  essentially  the  same. 

A  similar  calculation  can  be  made  for  the  isentropic  case,  curve  AC 
of  Fig.  5.  Then  elastic  constants  are  isentropic  instead  of  isothermal  and  T 
must  be  calculated  in  order  to  determine  entropy  delivered  to  the  reservoir. 

In  this  case  entropy  of  the  reservoir  increases  monotonical ly.  Temperature 
is  a  state  function  of  S,  V,  Ye  according  to  Eqs.  (46)  and  (48),  as 

dT  =  (3T/9S)  a  dS  +  (3T/3V)  dV  +  (9T/9Ye)  dYe  . 

V,Ye  S,Ye  S,V 

Set  dS  =  0  for  isentropic  compression  and  use  the  equalities  among  cross¬ 

derivatives  Of  E  to  obtain 

dT  =  -(3p/3S)  0  dV  +  (8V/3) ( 9x/9S)  dYe  .  (79) 

V,Ye  V,Ye 

It  will  now  be  found  that  T  goes  through  a  maximum  and  a  minimum  during  each 
cycle  and  returns  to  its  initial  value  at  the  end  of  a  cycle,  since  dT  is  an 
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+  .05 


-.05  -.05 


-.05  -05 

Figure  6  Strain  and  stress  variations  for  isothermal  cycle  of  linear 
elastic-plastic  material  with  constant  yield  stress.  (See 
Appendix  D.)  Abscissa  is  / | da  | 
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Wp  S(SAMPLE) 


S(RESERVOIR) 


Figure  7  Entropy  and  plastic  work  for  the  example  of  Figure  6 


exact  differential.  Since  the  dependence  of  observables  on  entropy  is  not 
measured  directly,  this  case  appears  to  be  of  lesser  interest  than  the  iso¬ 
thermal  case. 

Heat  conduction  effects  are  not  normally  significant  in  stress  wave 
propagation,  so  compression  is  inclined  to  be  adiabatic.  V,  ye,  S,  T  do  not 
repeat  cyclically  as  stress  is  cycled  because  plastic  work  is  retained  in  the 
sample  as  heat.  One  can  ask  for  the  fate  of  mechanical  and  thermal  variables 
as  p  increases  to  its  maximum  value  and  returns  to  zero.  It  may  be  of  further 

A 

interest  to  inquire  what  processes  are  required  to  return  it  from  this  final 
state  to  its  initial  state,  but  this  is  not  normally  of  interest,  for  example, 
in  a  shock  wave  problem. 

In  the  adiabatic  case,  both  entropy  and  temperature  are  changed  with 
each  increment  of  plastic  work.  The  sample  is  isolated  from  the  universe  so 
the  only  source  of  entropy  is  the  plastic  work  done  on  the  specimen.  Then 

dS  =  dWp/T  .  (80) 

Combining  this  with  Eqs.  (67)  to  (69)  yields 

dT  =  (T/CVy)  [dWp/T  -  (sp/3T)VjYe  dV 

+  (8V/3)  (3x/3T)v^Ye  dye]  (81) 

4.2  Plastic  Work  Partially  Converted  to  Internal  Strains. 

Thermodynamic  Potentials 

The  admission  that  plastic  work  changes  material  properties  has  rather 
profound  effects  on  thermodynamic  concepts  though  the  magnitude  of  energy 
involved  is  not  very  large;  probably  not  more  than  5  to  10%  of  the  plastic  work 
is  stored  in  this  way,  and  the  total  may  be  much  less.  The  most  directly 
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observed  effect  of  internal  energy  storage  is  work-hardening;  in  some  materials 
the  yield  stress  may  be  doubled  by  plastic  work.  Although  this  effect  results 
from  partition  of  the  plastic  work,  it  is  in  itself  not  a  thermodynamic  effect 
and  so  must  be  taken  as  an  empirical  result. 

In  order  to  treat  this  case  thermodynamically,  we  still  consider  that 
at  each  point  on  a  deformation  path  the  material  is  in  a  reproducible  thermodyna¬ 
mic  state,  but  that  the  state  depends  on  an  internal  variable  as  well  as  the 
external  variables  of  the  previous  section.  The  physical  and  mechanical  char¬ 
acteristics  of  plastically  deformed  material  are  commonly  ascribed  to  the 
creation  and  entanglement  of  dislocations,  so  dislocation  parameters  are  used 
as  internal  thermodynamic  variables.  Dillon  and  Kratchovil  (1970)  have  used 
dislocation  density  and  dislocation  entanglement  as  two  internal  variables. 
However  the  present  state  of  understanding  of  dislocations  hardly  seems  to 
justify  the  use  of  two  variables.  Here  the  dislocation  density  will  be  taken 
as  the  sole  internal  variable.  Let  this  be  denoted  by  N,  the  total  dislocation 

length  per  unit  volume  of  material.  Denote  the  energy  stored  per  unit  length 

o 

of  dislocation  by  pr  =  b  y/2  (Nabarro,  p.  694).  Then  changes  in  internal  energy 
are  reversible  in  the  variables  S,  V,  ye,  N: 

dE  =  TdS  -  pdV  +  (8Vt/3)  dye  +  rdN,  (82) 

implying  that 


E(S,  V,  ye,  N) 

(83) 

A(T,  V,  ye,  N). 

(84) 

These  potentials  being  given,  the  change  in  entropy  resulting  from  isothermal 

compression  can  be  calculated  as  before  (dT  =  0): 

dS  =  (9S/9V)  dV  +  (9S/9ye)  dTe 

T >y  »N  T,V,N 

+  (3S/9N)  e  dY6  • 

T,V,y 
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Using  the  compatibility  relations  derived  from  A  this  can  be  written 

dS  =  (ap/aT)  dV  -  (8V/3)(3t/3T)  dYe 

V,Ye,N  V,Ye,N 

+  (ar/aT)  dN  (85) 

V,ye,N 

4.3  Compatibility  Relations 

dS  is  an  exact  differential  since  it  is  derived  from  the  state  func¬ 
tion,  Eq.  (84).  This  implies  that  the  coefficients  in  Eq.  (85)  satisfy  the 
following  conditions: 


(a2p/9Taye)v>N  = 

- (8/3) (3[V( 9t/9T)  e  ]  /  3V) 

V,ye,N  ye,N 

(86a) 

(a2p/aNaT)  0  - 

v,  ye 

-(a2r/aVaT) 

Y  »N 

(86b) 

(a2r/aTaye)VjN  = 

(8/3)  {S[V(3t/3T)  e  ]  /  3N} 

\/  M  \/  6. 

(86c) 

V  >y  jN  V ,yc 

Assume  that  yield  always  occurs  at  sufficiently  small  shear  that  shear 
stress  and  strain  are  linearly  related, 


t  =  2pye  (87) 

ap/3ye  =  0.  (88) 

Assume  also  that  pressure  is  independent  of  elastic  shear  strain.  While  this 
is  not  necessarily  true,  it  is  highly  plausible: 

ap/3ye  h  0.  (89) 

From  Nabarro  (p.  694)  one  finds  that 

r  »  b2  Vy/2 .  (90) 
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The  implications  of  Eqs.  (86)  to  (90).  are  derived  in  Appendix  I.  They  are 


(9P/9T) 

Y6,V,N 

=  f2(T,V) 

(91) 

(9(yV)/3T) 

Y*.V 

=  MT) 

,N  1 

(92) 

P  =  f (T>V)  . 

+  g(v,N) 

(93) 

yV  =  F(T)  + 

G(V,N) 

(94) 

Equations  (91)  to  (94)  are  useful  in  constructing  equations  of  state  for  the 
classes  of  materials  described  here.  For  example,  Eq.  (93)  limits  the 
Helmholtz  potential  to  the  form 

A(T,V,Ye,N)  =  -  / f ( T , V ) dV  -  / g(V ,N)dV  +  4>(T,N,Ye)  (95) 

where,  f,  g  and  <f>  are  arbitrary  functions  of  the  variables  indicated. 

The  possibility  of  deriving  relations  like  those  of  Eqs.  (91)  to  (95) 
is  a  consequence  of  the  assumption  of  uniaxial  strain,  on  which  this  entire  cal¬ 
culation  is  based.  Under  conditions  of  uniaxial  strain  it  is  possible  for 
density  changes  and  plastic  strain  to  occur  simultaneously.  In  uniaxial  stress 
this  is  not  the  case:  thermodynami c  potentials  do  not  depend  on  both  V  and  y , 
so  compatibility  relations  of  the  kind  derived  above  do  not  exist.  It  does  not 
seem  likely  that  going  to  a  more  general  process  than  uniaxial  strain  would  pro¬ 
duce  additional  information,  except  in  anisotropic  materials. 

4.4  Gibbs  Relation  Including  Stored  Energy 

The  Gibbs  relation  for  entropy  production  in  an  arbitrary  process  can 
be  derived  as  before.  Equate  the  increase  in  internal  energy  of  unit  mass  to 
the  sum  of  work  done  and  heat  added;  this  gives  the  analogue  of  Eq.  (42): 
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TdS  +  dWe  +  rdN 


dW  +  dQ 


(96) 


and  the  Gibbs  relation  becomes 

TdS  =  dW  -  dWe  -  rdN  +  dQ 

=  dWp  -  rdN  +  dQ  .  (97) 

Equation  (97)  is  a  formalization  of  the  basic  assumption  that  energy  stored  in 
dislocations  can  be  recovered.  Plastic  work  may  be  required  to  move  disloca¬ 
tions,  but  this  is  directly  associated  with  plastic  strain  and  appears  as  heat. 
Since  the  energy  of  dislocations  derives  from  local  strain  fields,  internal 
strain  could  equally  well  be  used  as  a  parameter.  N  is  more  convenient  since 
extensive  work  exists  on  the  properties  of  dislocations  and  their  behavior. 

When  N  is  sufficiently  large,  dislocations  interact  and  their  energy  is  no 
longer  proportional  to  N.  Such  interactions  can  be  partially  accounted  for 
by  letting  r  be  a  function  of  N.  But  we  shall  ignore  them  here  and  intro¬ 
duce  self-annihilation  and  annealing  of  dislocations  to  prevent  N  from  growing 
indefinitely. 

If  the  compression  process  is  carried  out  at  a  finite  rate,  Eq.  (97) 
gives  the  rate  of  entropy  production: 

TdS/dt  =  dWp/dt  -  rdN/dt  +  dQ/dt  (98) 

In  the  event  that  either  dWp/dt  or  dN/dt  is  rate-dependent,  this  will  lead  to 
entropy  production  beyond  that  given  by  Eq.  (85). 

4.5  Equilibrium  Values  of  N;  An  Anomaly 

In  order  to  complete  the  thermodynamic  description,  a  procedure  is 
required  for  calculating  N.  Two  routes  are  available  for  obtaining  the  required 
information.  One  is  to  continue  the  formal  thermodynamic  description  of  the 
process  and  so  to  determine  the  bounds  on  equilibrium  values  set  by  the 
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compatibility  relations.  Such  information  can  be  augmented  by  statistical 
calculations  for  determination  of  entropy  and  lattice  calculations  for  deter¬ 
mining  such  parameters  as  r.  The  latter  has  already  been  noted  in  the  pre¬ 
ceding  paragraph.  The  second  is  to  apply  some  of  the  formalism  of  dislocation 
theory  and  expressions  for  multiplication  rates  to  relate,  for  example,  dislo¬ 
cation  density  to  strain  history.  Both  of  these  routes  are  considered  here. 

We  start  with  the  thermodynamic  development. 

The  great  question  about  the  thermodynamics  of  dislocations  is  whether 
equilibrium  ever  exists  in  the  sense  that  the  appropriate  potential  is  mini¬ 
mized  with  respect  to  the  various  parameters.  Here  we  assume  that  an  equilib¬ 
rium  density  N  exists  and  we  find  how  it  is  related  to  other  material  para¬ 
meters.  For  this  purpose  we  consider  Eq.  (82)  and  its  companion  for  dA: 

dE  =  TdS  -  pdV  +  (8V  /3)  dYe  +  rdN  (82) 

dA  =  -SdT  -  pdV  +  (8V  /3)  dYe  +  rdN.  (99) 

Because  Eqs.  (82)  and  (99)  are  exact  differentials,  the  cross  derivatives  of 
their  coefficients  must  be  equal.  Those  involving  N  are 


(5T/3N)  e  = 

s,v,Ye 

(9r/9S) 

v,y,n 

(100a) 

(9p/3N) 

-( 9r/9V) 

(100b) 

S,V,Ye 

S,Ye,N 

(16VYe/3)  (9y/9N) 

=  (9r/3Ye)s  v  N 

s,v,Ye  b,v,N 

(100c) 

(9S/3N)  = 

T,V,Ye 

-(9r/9T) 

V ,Ye,N 

- 

(b2V/2)  (9P/9T)  p 

(1 OOd) 

V,Ye,N 
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-(3p/aN)TjV^e  =  (3r/3V)y^e^ 

(16Vye/3)  (9y/3N)TjVjYe  =  (Sr/^e)T>V,N 


(lOOe) 

(lOOf) 


Here  the  approximations  of  Eqs.  (87),  (88)  and  (90)  have;  been  used. 

To  see  the  utility  of  Eqs.  (100),  consider  Eq.  (100a).  The  right 
hand  side  can  be  written 


(ar/aS)  =  (b^V/2)  ( 3y/aS) 


V,ye,N 


=  (b^V/2)  (3y/3T)  (3T/3S)  p 

V,ye,N  V,ye,N 


=  (b‘VT/2Cv  „)  (3y/3T) 

VyN  V,ye,N 


where 


CVyN  s  T<3S'9T>„  e 


V  ,y  ,N 


(101a) 


(101b) 


By  regarding  N  as  a  function  of  (T,  V,  y  ,  p),  we  have 


(3N/3T)  p  =  (3N/3T)  +  (aN/ap)  .  (3P/3T) 

S,V,ye  V,yS,p  T,V,y6  S,V,ye 


It  seems  quite  plausible  that  p  is  insensitive  to  N*,  then  the  second  term 
is  negligible  and 


(3N/3T) 


S,v,y* 


( aN/aT) 


V,Ye,P 


( 1 01  c ) 


Equations  (100a),  (101a)  and  (101c)  combine  to  give 


ON/a-O  e  -  2C  /  [b  vT( 3n/aT)  e  ] 
V,y  ,p  v,y  ,N 


(102) 
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This  is  a  curious  and  unexpected  result.  3y/3T  is  known  to  be  negative  and 
SyN  1S  Positive>  therefore  aN/aT  <  0.  But  one  expects  the  equilibrium 

number  of  dislocations  to  increase  with  T  for  the  following  reasons.  Energy 
is  required  to  form  dislocations,  so  that  an  assembly  which  includes  disloca¬ 
tions  contains  eigenstates  of  higher  energy  than  a  perfect  lattice.  Under 
equilibrium  conditions  these  higher  energy  states  are  unoccupied  at  low  tem¬ 
peratures  but  at  sufficiently  high  temperatures  they  may  be  occupied.  Just 
where  the  argument  leading  to  Eq.  (102)  fails  is  not  clear,  but  it  is  clear  that 
the  reason  for  the  failure  must  be  determined  before  other  relations  obtained 
from  the  compatibility  condition  can  be  safely  used. 

Changes  in  Cy^  with  N,  Ye  and  V  can  be  obtained  from  the  expression 
for  dS  when  dT  /  0.  Equation  (85)  is  then  replaced  by 

dS  =  (CVYN/T)dT  +  (ap/aT)dV  -  (8V/3)  (3t/3T)  dYe  -  (ar/3T)dN.  (103) 

Then 

3(cvyn/t)/3N  =  ‘  (104a) 

3(cVyN/t)/3Y6  =  -(8V/3)  32t/3T2  (104b) 

3(CVYN/T)/3V  =  O04c) 

If  Eqs.  (87)  and  (90)  are  assumed  true,  Eqs.  (104a)  and  (104b)  combine  to 
yield  the  result 

aCvyaN  -  ( 3b2/ 18)  9CVYN/a(Ye)2  =  0  .  (105) 

This  has  solutions  of  the  form 

CVYN  =  f[(N  +  6Yfi2/b2),  V,T],  (106) 
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which  shows  an  interesting  relation  in  the  dependence  of  on  N  and  ye. 
It  does  not  appear  to  be  susceptible  to  test  with  present  information. 

The  physical  consequences  of  relations  like  those  of  Eqs.  (100)  and 
(104)  have  been  illustrated  but  not  exhausted.  For  the  present,  however,  we 
forsake  further  development  of  these  relations. 


48 


V.  MECHANICAL  CALCULATIONS  OF  DISLOCATION  DENSITY,  N 
5.1  Formalism 

As  a  consequence  of  plastic  work,  dislocation  densities  may  be  generated 
which  are  much  greater  than  equilibrium  values  calculated  in  the  previous  section. 
These  may  represent  substantial  quantities  of  stored  energy,  and  at  sufficiently 
high  temperatures  this  should  be  released  as  heat  with  a  consequent  reduction  of 
N.  The  thermodynamic  effects  of  non-equilibrium  changes  of  N  are  derived  from 
the  Gibbs  relation  in  the  form  of  Eq.  (100).  This  indicates  that  dN/dt  is  the 
quantity  of  prime  interest. 

A  form  frequently  used  for  dN/dt  is  quoted  by  Webster  (1966): 

dN/dt  =  aN  -  bN2  +  SQ.  (107) 

The  first  term  represents  creation  of  dislocations  by  regenerative  processes, 
the  second  represents  self-annihilation,  and  the  third  represents  spontaneous 
nucleation  of  dislocations  under  the  influence  of  externally  applied  forces. 
Equation  (107)  can  be  integrated  directly  to  yield 

N  =  [A  -  Be  exp(-Kt)]  /  [1  +  e  exp(-<t)]  (108) 

where  k  =  (a2  +  4b  Sq)^2 

A  =  (a/2b)  (1  +  K/a) 

B  =  -(a/2b)  (1  -  K/a) 

0  =  (A  -  N.)/(B  +  N.) 

N.j  =  dislocation  density  at  t  =  0 

For  SQ  e  o, 

N  =  a  N.  exp(at)  /  [a  +  N^b  (exp  (at)  -  1)  ].  (109) 
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The  Gilman  relation  for  multiple  cross  glide  has 

aN  =  mbv  N  =  N.  f  yP 

where  yP  is  plastic  strain  rate.  Then,  for  small  multiplication,  changes  in 
N  are  proportional  to  plastic  strain.  This  relation  has  been  widely  used 
in  discussion  of  dynamic  failure  (Asay,  et  al.  1972). 

For  large  t  and  Sq  =  0,  Eq.  (109)  goes  to 

N  =  N$  =  a/b  (110) 

This  is  the  saturation  value  of  N;  it  is  unrelated  to  the  thermal  equilibrium 
value  discussed  in  the  previous  section.  From  Eq.  (109),  1/a  =  t0  is  a 
characteristic  time  for  equilibration.  Then  Eq.  (110)  can  be  rewritten  in 
the  form 

1/N  -  1/N$  =  (1/N.  -  1/N$)  exp(-t/te);  (111) 

t  is  assumed  to  be  the  order  of  a  few  multiplication  times  and  is  discussed 
e 


in  Appendix 

H: 

V 

= 

ai/v.  -  (a/n)/  £n  [1  -  exp(-U/kT)] 

'O 

(112) 

where  6 

= 

distance  between  imperfections  which  cause 

multiplication 

n 

= 

a  characteristic  frequency  associated  with 

pinned  dislocations 

U 

= 

Uq  (1  -  o/oq)  is  an  activation  energy 

V4 

= 

shear  wave  velocity 

a 

= 

a  number  the  order  of  ten 

a 

= 

effective  stress  driving  the  dislocation 

a0 

critical  or  threshold  stress 

. 

Equation  (112)  is  to  be  applied  if  o.<  cjq.  For  a  >  Oq, 

t0  =  ai/v^.  (113) 
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If  U/kT  'v-  10,  Eq.  (112)  gives  t£  10"^  seconds,  which  is  very  long  on  the 

scale  of  shock  experiments.  On  the  other  hand,  Eq.  (113)  is  the  order  of 
-12 

10  for  a.100  pinning  points  per  million  lattice  sites,  a  time  which  is 
completely  negligible. 


5.2  Applications  to  Shock  Processes 

If  we  assume  the  Gilman  relation  for  multiplication  and  ignore 
self-annihilation,  we  have 


or 


dN/dt  *  mbvN  =  my*3 
dN  =  mdyP. 


Each  increment  in  plastic  work  is  also  proportional  to  dyP: 

dWp  =  (8VT/3)dYp. 

The  increment  in  stored  energy  is  rdN  and  its  ratio  to  dWp  is 

rdN/dWp  =  3mr/8 Vt  =  3b2ym/16x  . 

For  m  =  1011,  b  =  3  x  10-8,  y  =  2.7  x  1011  and  r  =  5  x  109  dynes/cm2, 
this  ratio  is 

rdN/dW  *  .9  x  10"3  =  .09% 

This  value  is  to  be  compared  with  the  oft-quoted  value  of  "5  to  10%."  This  can 
be  achieved  by  making  t  very  much  smaller,  as  in  a  static  test,  or  by  increasing 
m.  A  value  of  0.5  kilobars  in  a  static  test  on  soft  aluminum  is  not  out  of 
reason,  in  which  case  the  ratio  of  stored  energy  to  plastic  work  goes  to  0.9%. 
Moreover  m  is  not  precisely  known  and  may  be  increased  several -fold  without 
seriously  violating  any  experiments.  In  fact  some  of  the  dynamic  failure  data 
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12 

for  aluminum  suggest  that  m  may  be  as  large  as  2.5  x  10  in  shock  experi¬ 
ments.  Then  taking  x  =  5  kb  again,  . 

rdN/dWp  =  2.3% 

which  is  approaching  the  usually  quoted  range  for  stored  energy.  The  possibility 
that  plastic  stored  energy  in  shock  compression  is  less  than  static  values  must 
be  admitted. 

For  comparison,  a  similar  calculation  for  the  case  of  uniaxial  stress 
leads  to  the  ratio 

rdN/dWp  =  mr/VY  =  mb2p/2Y 

where  Y  =  2x  is  the  yield  stress  in  simple  tension.  This  does  not  differ 
significantly  from  the  value  for  uniaxial  strain. 

The  total  plastic  work  generated  when  aluminum  is  shocked  to  lOOkb  is 

O 

approximately  Wp  =  .15  kb  cc/g  or  1.5  x  10  ergs/g.  Equating  5%  of  this  to 

_5 

rN  where  r  =  4.5  x  10  gives 

N  =  .05  x  1.5  x  108/4.5  x  10"5  =  1.6  x  1011. 

This  is  small  for  a  saturation  value  of  N  and  it  suggests  that  in  the  small 
times  available  under  shock  conditions,  saturation  of  dislocation  density  may 
not  occur. 
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VI.  DISCUSSION  AND  CONCLUSIONS 


A  simple  and  approximate  model  for  treating  the  thermodynamics  of 
elastic-plastic  solids  has  been  initiated  here  but  not  completed.  Its  appli¬ 
cation  to  conditions  of  unixial  strain  shows  that  the  additional  thermodynamic 
degree  of  freedom  compared  to  uniaxial  stress  gives  some  potentially  interesting 
new  compatibility  relations  which  can  be  checked  against  various  experimental 
data.  It  is  also  suggested  that  shock  experiments  which  have  uniaxial  strain 
geometry  may  be  useful  in  evaluating  dislocation  concepts  and  those  of 
el asti c-pl as ti c  thermodynami cs . 

The  point  of  view  espoused  here  is  a  novel  one.  It  is  common  to  speak 
of  an  elastic-plastic  system  as  being  "essentially  irreversible"  or  "immersed 
in  a  sea  of  irreversibility."  The  view  expressed  here,  but  not  completely 
developed,  is  that  the  elastic-plastic  state  is  completely  reversible  in  the 
variables  which  are  natural  to  the  state  and  that  the  apparent  irreversibility 
associated  with  inability  to  retrace  a  o-e  state  is  an  awkward  consequence 
of  geometry  and  is  not  essentially  related  to  the  question  of  reversibility  of 
the  process.  Having  stated  such  a  hypothesis,  one  must  then  make  some  decision 
about  the  role  of  stored  plastic  work.  For  many  purposes  the  problem  is 
indifferent  to  this  decision  because  the  quantity  of  stored  energy  is  small. 

But  when  the  stored  energy  is  related  to  dislocation  density,  the  thermodynamic 
formalism  suggests  that  this  may  provide  some  new  avenues  for  investigation  of 
the  role  of  dislocations  in  plastic  deformation. 
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APPENDIX  A 


Addition  of  Elastic  and  Plastic  Strain  Increments 

In  Fig.  A1  are  shown  two  points  of  a  continuum,  0  and  P,  separated  by  a 
small  distance  s.  Application  of  an  external  stress  forces  0  to  A  and  P  to  B. 
When  the  stress  is  removed,  A  relaxes  to  Q  and  B  to  T.  We  have  from  the  figure 


fflt  V 

-> 

U, 

— y 

PB 

=  u 

+  du 

du 

-> 

=  s 

•  vu 

(Al) 

AQ  = 

U’, 

BT 

=  u1 

1  +  du' 

du' 

=  (< 

;  +  du)  •  vu1 

(A2) 

OQ  = 

u". 

— * 

PT 

=  0" 

+  du" 

du" 

-y 

=  s 

•  vu" 

(A3) 

In  writing  Eqs.  (Al)  -  (A3)  the  vector  change  du,  etc  has  been  expanded  in  a 
Taylor  series  in  s  or  s  +  du  and  all  but  the  lowest  order  term  discarded. 

du"  is  identified  as  the  plastic  displacement,  -du'  the  elastic  dis¬ 
placement  and  du  as  the  total  displacement.  Then 

-du 1  +  du"  =  s  •  (-vu 1  +  vu")  -  du  •  vu' 

=  s  •  (-vu 1  +  vu"  -  vu  •  vu ' )  (A4) 

The  vector  gradient,  vu,  can  be  decomposed  into  a  strain  tensor  4  and  a  rotation 
tensor  6: 

VU  =  £  +  £ 

where  4  =  (1/2)  (vu  +  vu) 

£  =  (1/2)  (vu  -  vu) 

and  vu  is  the  transpose  of  vu. 

Then  Eq.  (A4)  becomes 

-du'  +  du"  -  s  •  (-4'  +  4"  -  4'  •  4"  -  4'  *  4"  +  rotation  terms). 

Then  if  e.' ,  4",  0 ' ,  4"  are  infinitesimal,  the  2nd  order  terms  can  be  neglected  so 
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that 


du  =  -du'  +  du"  =  s  •  +  4"  +  rotation  terms)  =  s  •  (£  +  g) 

so  that 

e_  =  4"  =  +  (A5) 

where  ef  has  been  set  equal  to  -4'.  In  a  finite  deformation  process  the  relation 
can  be  applied  to  each  increment  in  strain  so  that 

de_  =  dg  +  dej3 .  (A6) 


Figure  A1  Arbitrary  displacements  in 
a  continuum 
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APPENDIX  B 


Stresses  in  Elastic-Plastic  Deformation 

The  conventional  assumption  in  elastic-plastic  theory  is  that  stress  is 
supported  by  the  elastic  strain  and  there  exists  no  contribution  from  the  plastic 
strain.  This  appears  plausible  if  one  examines  the  microscopic  behavior  of  a 
plastically  deforming  material  from  the  point  of  view  of  dislocation  theory. 

Plastic  deformation  is  synonymous  with  motion  and  generation  of  disloca¬ 
tions.  In  most  materials  dislocations  probably  move  at  all  stress  levels,  so 
there  is  no  such  thing  as  a  yield  point  and  there  is  no  truly  elastic  behavior. 
But  in  practice  the  yield  point  concept  is  useful  and  the  yield  point  itself 
may  be  taken  to  be  the  point  at  which  large  numbers  of  dislocations  are  set  in 
motion.  The  motion  of  dislocations  is  inhibited  by  the  existence  of  energy 
barriers  which  must  be  overcome  by  combination  of  stored  elastic  energy  and 
thermal  fluctuations  of  atoms.  So,  for  example,  one  may  consider  a  strained, 
work  hardening  solid  as  one  in  which  many  dislocations  exist  but  are  momentarily 
immobilized  or  "pinned"  by  energy  barriers.  If  applied  stress  is  held  constant, 
dislocations  may  occasionally  overcome  a  barrier  and  move  to  the  next  barrier, 
thus  contributing  to  the  plastic  deformation  by  creep. 

If  the  external  stress  is  increased,  local  strain  energy  is  increased 
around  the  pinning  points,  more  dislocations  are  moved  through  pinning  points 
and  pass  on  until  they  are  pinned  again.  If  it  is  assumed  that  dislocations 
move  freely  between  pinning  points,  the  plastic  strain  which  results  from  their 
motion  requires  no  part  of  the  applied  stress;  i.e.  the  applied  stress  is 
supported  entirely  by  the  elastic  part  of  the  strain.  It  does  follow,  however, 
that  a  relaxation  process  may  exist.  For  example,  an  increment  in  stress 
increases  local  strain  energy  around  pinning  points  and  increases  probability 
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that  dislocations  will  break  free  of  their  pinning  points  and  move  on  until 
pinned  again.  However,  the  probability  of  a  dislocation  breaking  free  is  time 
dependent  because  it  depends  on  thermal  fluctuations  as  well  as  strain  energy. 

This  means  that  an  increment  in  stress  will  produce  an  immediate  elastic  strain, 
and  that,  as  time  passes,  dislocations  escape  their  pins,  move,  and  are  repinned, 
so  that  a  plastic  strain  develops.  Or  the  original  increment  in  strain  is 
reduced  as  the  pinned  dislocations  escape  from  their  barriers  and  the  elastic 
strain  is  reduced.  This  concept  can  be  expressed  quantitatively  in  the  follow- 
i ng  way 

da  =  pa2  de6  (B1 ) 

=  pa2(de  -  deP) 

=  pa2de  -  a2(deP/dt)dt 

or  da/dt  =  pa2de/dt  -  Pa2  deP/dt. 


In  terms  of  stress  and  strain  deviators  this  translates  to 


dS./dt  =  2y  de./dt  -  2yde.p/dt 

J  J  J 


=  2y  de./dt 
J 


F(S,e) 


(B2) 


where  F  is  a  relaxation  function. 

Equations  (Bl)  and  (B2)  are  both  based  on  the  assumption  that  dislocations 
move  between  pinning  points  without  drag.  If  this  is  not  true,  then  a  viscous 
contribution  to  the  stress  appears: 

da  =  pa  +  n  de 


or  dS./dt  =  2y  de./dt  -  2y  de.p/dt  +  2nde.p/dt 

J  J  J  J 


(B3) 


where  e .p  =  de.p/dt. 

J  J 


60 


APPENDIX  C 


Plastic  Work  and  Plastic  Strain 
From  Eq.  (24b) 

dexP  =  dEx  •  dexe  <cl> 

Eq.  (24a)  d6  =  de  +  2de  =  de  e  +  2de  e  +  (de  P  +  2de P) 

Ay  a  y  x  y 

=  dexe  +  2dEye  +  0 

=  dex  by  Eq.  (28a). 


V  ■ 

dex/2  -  dexe/2  =  dexP/2 

d°x  = 

Ado  +  2y  de  e  by  Eq.  (27) 

A 

doy  = 

Ado  +  2y  d ey  ^  by  Eq.  (27) 

d(ax-ay)  = 

2y(dexe  -  deye)  =  2y(dex6  -  dexP/2) 

(02) 

dx  = 

[d(ax-oy)]/2  =  y  dexe  -ydexP/2 

dExe  ' 

dx/y  +  dexP/2 

(C3 ) 

dexP  * 

de  -  dx/y  -  de  P/2  from  Eqs.  (Cl)  and  (C3) 

A  A 

dExP  ■ 

(2/3)  (dex  -  dx/y) 

(C4) 

Combine  this  with  Eq.  (36)  to  obtain 

dW  =  t  Vt  •  T  de  P  =  2Vxde  P  (C5) 

P  W  dm  X  X 
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From  Eq.  (25)  dWp  =  V(SxdexP  +  2Sy  deyp) 

=  jVt  (dexp  -  deyp)  from  Eq.  (30)  and  (31)  (C6) 

dex  =  dex  -  de/3  =  dexe  +  dexP  -  de/3 

dey  =  dey  -  de/3  =  deye  +  deyP  -  de/3 

dex  '  dey  =  dex  -  dey 

(dexe  -  dey6)  +  (dexP  ■  deyP)  =  (dexe  "  deyG)  +  (dexP  "  deyP)  (C7) 

dexe  -  deye  =  dT/y  by  Eq.  (32)  ) 

1  (C8) 

dex6  "  dey6  =  dx/^  by  Eq.  (C2)  ) 

Therefore  Eq.  (C7)  becomes 

de/-de/  -  Ce/.dc/ 

since  de  p  +  2de  p  =  0. 

*  y 

Substituting  this  into  Eq.  (C6)  yields 

dWp  =  y  Vt  •  dexP  =  2VrdexP 
in  agreement  with  Eq.  (C5). 

Note  that  the  following  argument  is  faulty: 

dex  ■  dEx  -  d9/3  “  I  d£x 

dey  *  dcy  '  de/3  “  7  dEx 
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dex  -  dey  =  dex 


(dexe  -  deye)  +  (dexp  -  deyp)  =  de/  +  d£)(e 
•*  dexp  -  deyp  =  dexp 
and  dexe  -  deye  =  dexe 

This  result  contradicts  Eq.  (C8)  and  is  not  true  because  the  non-zero  values- of 
de  e  and  de  p  have  been  ignored  in  obtaining  the  final  relation. 

y  y 
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APPENDIX  D 


Stress  Cycles  for  a  Linear  Elastic-Plastic  Material  with  Constant  Yield  Stress. 
All  Plastic  Work  is  Converted  to  Heat. 


T  =  300°K 

A  =  constant  =  582  kb 

y  =  constant  =  274  kb 

Y  =  constant  =  10  kb 


e 


£ 


z 


-V ( 9o  / 3T)  p  =  16.8  x  10~3  kbcc/g° 
x  £xex 


VQ  =  .37  cm3/g 

Cv  =  8.4  x  10"3  kbcc/9° 

r  =  2 


Initial  Conditions  at  0  (Figs.  3  and  4) 


Wp  =  0,  S  =  0;  entropy  of  heat  reservoir,  =  0 


Along  0A  (Figs.  3  &  4) 
d£x  =  dE^6  =  dax/(A  +  2y ) 


day  =  Adc?x/(A  +  2y) 

d£  P  =  de  P  =  dW„  =  de  ,e  =  0 

x  y  P  y 

dWp  =  0 

dS  (reservoir)  =  -dS  (sample) 


Along  AB 

ax  -  ay  =  -Y;  d(ax  -  ay)  =  -dY  =  2y(d£xe  -  d£ye)  =  0 
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d*x 

■  dExe  + 

2deye  =  3dexe 

dox 

=  3xde  e  + 

A 

Zyckx6 

dex6 

=  da/ (3A 

A 

+  2r) 

dex 

■  "S*  ■ 

3dax/(3A  +  2y ) 

d^xP 

■  dex  •  d£xe  -  V 

dwp 

=  2Vtde  P  ■ 

X 

■  *VVdExP 

V  = 

v0  exp(ex) 

day 

=  da 

X 

deyP 

=  -se 

dS  (reservoir)  =  -dS  (sample)  +  dW^/T 
Along  BC 

dex  =  dexe  =  dax/(x  +  2y) 

day  =  Xdax/(A  +  2y ) 
dExP  =  deyP  *  dey6  *  0 

dWp  •  0 

dS  (reservoir)  =  -dS  (sample) 
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Aong 

_CF 

°x  ‘ 

°y  =Y 

dey6 

*  d*x" 

d£x 

■  V  ■ 

3dax/(3x  +  2y) 

d£xP 

■  2dV3 

deyP 

*  -dox/(3x 

+  2y ) 

day 

=  dox 

dwp 

-  VV  dexP 

dS  (reservoir)  =  -dS  (sample)  +  dWp/T 


Alone 

l  FO 

dEx 

=  de/  = 

dox/(X  +  2 y) 

doy 

-  Xdax/(X 

+  2y) 

dExP 

=  de  P  • 

y 

■  se  =  0 

dW 

p 

=  0 

dS  (reservoir)  =  -dS  (sample) 
At  0 
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Compressive  stress  at  B:  assume  c„  =  -40.6  kb 

A  L 


a  »  a 


St  A:  °x  "  °y  =  n  -  x/ ( X  +  2y)]ax  =  -Y 

a„A  =  -U  +  2y)Y/2y  =  -20.6  kb;  a.A  =  -10.6  kb 
*  _ _ _ 


at  B: 


a  B  =  -40.6  kb  (assumed),  oB  =  -30.6  kb 

*  JL . . . . 


at  C:  a  C  -  a  C  =  Y 

a  y 


C  C  /  B  Bv  ,  C  B%  /  C  Bs 

aX  ‘  °y  "  (ax  ■  °y  )  =  (ax  “  ax  )  ■  (°y  "  °y  ) 


(oxC  -  oxB)  [1  -  x/(x  +  2y)]  =  2Y 


axC  -  0XB  =  V(X  +  2y)/y  =  41.2  kb 


av°  =  0.6  kb 

A 


■9.4  kb 


F  F 

at  F:  °x  "  °y  =  Y  whether  the  point  F  is  reached  by  expanding  along  the 
elastic  isotherm,  OF,  or  along  the  plastic  curve,  CF.  In  the  former  case, 


axF  =  -a/  =  20.6  kb.  Then  a  F  =  10.6  kb 

■  - . .  . ■  ■  .  — mL  .  „  ■—  


at  0:  Both  a x  and  return  to  zero  since  both  have  their  elastic  values  at  F 


eXA  =  °XA/(X  +  2u)  =  --01823 
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ewB  =  e.A  +  3(awB  -  a„A)/(3A  +  2y)  =  -.01823  -  .02616  =  -.04439 

a  A  A  A  *'*”"*■ 

e.,C  =  e.,B  +  (awC  -  a„B)/(x  +  2y)  =  -.04437  +  .03646  =  -.00793 

A  AAA  *w^™*m*” 

e..F  =  e..C  +  3{a.,F  -  o.,C)/(3A  +  2y)  =  -.00793  +  .026156  =  +.01823 

A  A  A  A 

£  0  =  Q 

X 

e 

fx_ 

’  ex6(A)  =  e/  =  -.01823 

e  e(B)  =  e  e(A)  +  (oB  -  a  A)/( 3A  +  2y)  =  -.01823  -  .00871  =  -.02694 

A  A  A  A 

eve(0  =  e  e(B)  +  (a  C  -  aB)/(A  +  2y)  =  -.02694  +  .03646  =  .00952 

c  e(F)  =  e  e(C)  +  (a/  -  avC)/(3A  +  2y)  =  .00952  +  .00871  =  .01823 

A  A  A  A 

£xe(0)  =  0 

fZ 

exP(A)  =  0 

exP(B)  =  2(axB  “  °xA)/{3X  +  2p)  =  --017428 

SP(0  =  .>) 

exP(F)  =  £XP(C)  +  2(  0XF  -  oxC)/(3l  *  2p)  =  0  =  exP(0) 

e 

ey6(A)  =0 

eye(B)  =  (axB  "  axA)/(3A  +  2v)  =  --008714 
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sye(C)  =  eye(B) 

ey6(F)  =  ey6(C)  +  (axF  "  aXC^^3X  +  2Vl)  =  0  =  ey6(0) 


Ye( A)  =  -.01823/2  =  -.00911 

ye(B)  =  -.00911 

ye(C)  =  +.00911 

ye(F)  =  .00911 

y6(0)  =  0.0 

Wp(A)  =0.0 

Wp(B)  =  -(2YVo/3)  [exp(£xB)  -  expUxA)]  =  .0625 

Wp(C)  =  Wp(B) 

Wp(F)  =  Wp(B)  +  (2YVq/3)  [exp(exF)  -  exp(e>(C)]  =  .0625  +  .0648  = 


.1273 
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Computation  of  Entropy 


From  Eqs.  (66)  -  (68)  changes  in  entropy  of  the  sample  are  given  by 

dS  =  rCv  dex  -  (16V/3)  (3y/ar)  e  Ye  dYe  . 

y  ,V 

From  Eq.  (35),  dYe  =  dx/2y. 

Since  x  is  assumed  constant  during  plastic  flow,  dYe  =  0  on  AB  and  CF.  Also, 

0  0.0  0 

dY  =  (de  -  de  )/2.  In  the  absence  of  plastic  flow,  dev  =  de  and 
x  y  xx 

de  e  =  0,  so  dYe  =  de  /2  on  0A,  BC  and  FO. 
y  x 

From  Appendix  J,  for  aluminum 

-(3y/3T)  _  *  .14  kb/°C  at  300°K. 

v,Ye 

From  Eq.  (75a),  V(9y/9T)  is  constant.  For  V  =  VQ  =  .37cc/g,  we  have 

- ( 1 6V/3) ( 9y/3T)  -  .27  kb  cc/g°C. 

Y  » V 

For  r  =  2,  Cy  =  8.4  x  10~3  kb  cc/g°C,  we  get 

dS  =  (16.8  x  10"3  dex  +  .27  YedYe)  kb  cc/g°C. 

From  0  to  A,  dYe  =  de  /2  and  Ye  =  e  /2,  so  at  A 

X  X 

S(A)  =  16.8  x  10"3  £  (A)  +  .03375  e..2(A)  =  -2.950  x  10"4. 

X  X 

On  AB ,  dYe  =  0,  so 

S(B)  =  S(A)  +  16.8  x  10“3  [e  (B)  -  e  (A)]  =  -7.345  x  10"4 

X  X 

From  B  to  C  dYe  =  de/2  but  Ye  /  e  / 2,  so  it  is  simpler  to  retain  Ye. 

X  X 
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Thus 


dS  =  16.8  x  10"3  de  +  .135  d(Ye)2 
S(C)  =  S(B)  +  16.8  x  10-3  [e  (C)  -  ev(B)] 

A  X 

+  .135[ye  (C)  -  ye2(B)]  =  -1.220  x  10~4 

From  C  to  F  dYe  =  0,  so 

S(F)  =  S(c)  +  16.8  x  10'3  [ex(F)  -  ex(C)]  =  3.175  x  10~4 
Finally  at  0 

S(0)  =  S ( F)  +  16.8  x  10’3  [e  (0)  -  e  (F)] 

A  X 

+  .135  [(Ye)2  (0)  -  (Ye)2  (F)]  =  4.68  x  10~9  *  0.0 


Entropy  of  the  Reservoir 
Sr(A)  =  -S{A)  =  2.950  x  10"4 

Sr(B)  =  Sr(A)  -  (SB  -  SA)  +  Wp(B)/300  =  9.429  x  10~4 

SR(C)  =  Sr(B)  -  (Sc  -  SB)  =  3.303  x  10~4 

Sr(F)  =  SR(C)  -  [Sp  -  Sc]  +  [Wp( F)  -  Wp(C)]/2  =  1.068  x  10~4 

SR(0)  =  SrF  -(S0  -  Sp)  =  4.243  x  10~4 
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TABLE  D-l.  Values  of  mechanical  and  thermal  variables  at  various  points  on 
the  cycle  of  Fig.  4.  Calculated  from  Program  "CYCLE."  Columns 
A,  B,  etc.  correspond  to  points  A,  B,  etc.  in  Fig.  4. 

First  entry  in  each  box:  Isothermal  cycle,  T  =  300°K 

Third  entry  in  each  box:  Isentropic  cycle,  S  =  0 

Fourth  entry  in  each  box:  Adiabatic  cycle 

The  entropy  in  parenthesis  is  obtained  from  exact  integration  as  described  in 


the 

first  part  of 

this  Appendix. 

A 

B 

C 

F 

0 

ov 

-20.6206 

-40.6000 

.642331 

20.6192 

0.0 

X 

(-20.6) 

(-40.6) 

(0.6) 

(20.6) 

(0.0) 

-20.5921 

-40.6000 

.461812 

20.6454 

0.0 

11 

II 

.380745 

20.3785 

it 

-10.6206 

-30.6001 

-9.35767 

10.6193 

-.000496 

y 

(-10.6) 

(-30.6) 

(-9.4) 

(10.6) 

(0) 

-10.5921 

-30.6001 

-9.53819 

10.6455 

-.000089 

II 

II 

-9.61926 

10.3786 

-.096066 

p 

13.9540 

33.9334 

6.02433 

-13.9526 

.000330 

13.9254 

33.9334 

6.20485 

-13.9788 

.000059 

II 

II 

6.28592 

-13.7119 

.064044 

V 

.363308 

.353937 

.367091 

.376806 

.369991 

.363308 

.353852 

.367002 

.376806 

.369991 

II 

.353837 

.366987 

.376736 

.369987 

£ 

-.0182483 

-.0443763 

-.00787783 

.0182470 

-4  x  10"8 

X 

(-.01823) 

(-.04439) 

(-.00793) 

(.01823) 

(0.0)  o 

1  x  10-% 

-.0182482 

-.0446184 

-.00812133 

.0182478 

II 

-.0446597 

-.00816250 

.0180612 

-1.1  x  10  s 

e 

Y 

-.00912418 

-.00912418 

.00912448 

.00912448 

.9  x  10"6 

(-.009113) 

(-.009113) 

(.009117) 

(.009117) 

(o.o)  7 

-.00912413 

-.00912413 

.00912403 

.00912403 

1  x  10  ' 

II 

II 

.00912406 

.00912406 

t 

l 

o 

X 

OO 

CO 

e  p 

0.0 

-.0174187 

-.0174187 

-2  x  10'6 

-2  x  10"6 

X 

(0.0) 

(-.017428) 

(-.017428) 

(0.0)  7 

(0.0)  7 

0.0 

-.0175801 

-.0175801 

-6  x  10-/ 

-6  x  10_/  - 

II 

-.0176076 

-.0176076 

-.000125060 

-1.3  x  10"4 

e 

ev 

-.0182483 

-.0269576 

.00954075 

.0182490 

2  x  10"6 

X 

(-.01823) 

(-.02694) 

(.00952) 

(.01823) 

(o.o}7 

-.0182482 

-.0270382 

.00945874 

.0182484 

6  x  10  ‘  . 

II 

-.0270520 

.00944509 

.0181863 

1.1  x  10" * 
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TABLE  D-l  —  Continued 


A 

B 

C 

F 

0 

0.0 

-.00870935 

-.00870935 

-1  x  10"6 

-1  x  10"6 

(0.0) 

(-.008714) 

(-.008714) 

(o.o)  7 

(0.0)  ' 

0.0 

-.00879006 

-.00879006 

-3  x  10_/  , 

-3  x  10  '  , 

II 

- . 00880379 

-.00880379 

-6.25  x  10'* 

-6.2  x  10"' 

-.000306570 

-.000745514 

-.000132340 

.000306558 

7  x  10“9q 

(-.000295) 

(-.000735) 

(-.000122) 

(.000317) 

(5  x  10"y) 

0.0 

0.0 

0.0 

0.0 

0.0 

ll 

+.000195374 

.000195374 

.000406573 

.000406573 

0.0 

.0624640 

.0624640 

.127240 

.127240 

(0.0) 

(.0625) 

(.0625) 

(.1273) 

(.1273) 

0.0 

.0630353 

.0630353 

.128409 

.128409 

II 

.0631326 

.0631326 

.128138 

.128138 

300 

300 

300 

300 

300 

(300) 

(300) 

(300) 

(300) 

(300) 

311.145 

327.989 

304.891 

289.223 

299.966 

II 

335.734 

312.090 

303.679 

314.849 
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Program  CYCLE 


In  the  general  case  where  X,  y,  Y  may  be  functions  of  density  and 
temperature,  numerical  integration  of  the  preceding  equations  may  be  required. 
This  can  be  accomplished  with  the  FORTRAN  program  described  here.  With  it  one 
can  obtain  values  of  the  variables  listed  above  for  isothermal,  isentropic  or 
adiabatic  cycles.  The  last  of  these  is  a  cycle  only  in  the  sense  that  stress 
returns  to  its  original  value.  The  other  variables  do  not  because  plastic  work 
is  retained  in  the  sample. 

Results  of  integrations  for  conditions  described  in  the  first  part  of 
this  Appendix  are  listed  in  Table  D-l  for  the  vertices  ABCFO.  For  isothermal 
and  isentropic  cases  the  cyclic  variables  do  not  quite  return  to  zero.  These 
integrations  are  accomplished  by  a  direct,  approximately  second-order  integra¬ 
tion  scheme.  The  listed  values  were  obtained  with  da  =  -.406  kb;  reducing 

A 

it  to  a  half  or  a  third  of  this  value  did  not  significantly  improve  the  closing 
values.  For  greater  accuracy  a  better  integration  scheme  would  be  required. 

The  following  assumptions  are  made  about  the  equation  of  state: 
bulk  modulus: 

K(T,V)  =  Kq  +  bp  +  6(T-T0)  +  gp(T-T0) 

where  K  =  isothermal  bulk  modulus  at  T  =  T  ,  p  =  0 
0  0 

b  =  polytropic  coefficient,  i.e.  p(TQ,V)  =  (KQ/b) [(VQ/V)b  -1] 

6  =  6  -  a  Kq  b 

a  =  thermal  expansion  coefficient,  (!/V)(aV/aT)p 
3  =  Cv(rQ  +  r1)/VQ 

Cy  =  specific  heat  at  constant  volume 
VQ  =  specific  volume  at  T  =  TQ,  p  =  0 

r  =  Gruneisen  parameter  =  rQ  +  r-j(l-  V//Vo^ 

9  =  (3/3Kq)  (4-1 /b)  -  ab2 
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w  -  p0  +  b,  (T/V  -  T0/V0) 

Tq  =  const,  is  a  reference  temperature 

For  aluminum  these  constants  are  approximately 

ro  ”  2.14,  r1  =  -3.2,  Cy  =  8.4  x  10  3  kb  cc/g°K,  VQ  =  .37  cc/g 

a  =  .75  x  10'Vk,  Kq  =  A  +  2y/3  =  764.67  kb 

A  =  582  kb,  yQ  =  274  kb,  Tq  =  300°K 

6  =  -.024065  kb/°K,  6  =  -.26895  kb/°K 

g  =  -1.40697  x  10" 3/°K,  b  =  4.27 

Provision  is  also  made  for  introducing  work-hardening  in  an  elementary 
way.  The  yield  strength  can  vary  with  plastic  work  according  to  the  relation 

Y  =  Y0  +  A[1  -  exp(Wp/W0)] 

where  YQ,  A  and  WQ  are  constants. 
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Comments  on  the  Program  Listing 


IC  is  the  integer  index  which  determines  whether  the  computation  will  be  made 

under  isothermal,  isentropic  or  adiabatic  conditions. 

R  R 

N  The  stress  increment  used  in  integration  is  a  /N  where  -a  is  the 

A  /\ 

maximum  compressive  stress  to  be  reached. 

MM  This  is  the  number  of  cycles  to  be  calculated.  For  this  run  two  cycles 
are  calculated.  IC  =  2  for  the  first  and  IC  =  3  for  the  second.  IC  is 
incremented  by  card  number  0184  (left  hand  column). 

m  = 


K0  =  K 


B  =  b 


GAMO  =  r 


GAM!  =  r1 


C  V  =  C, 


BETA  =  3 
DELTA  =  6 
SGAM  =  g 
AMU  =  b1 

Y0  =  initial  value  of  the  yield  stress 
A  =  A 


W0  =  WQ 

D 

SIGXB  =o  =  -(compressive  stress  at  which  stress  cycle  is  reversed) 

A 

MCOUNT  =  integer  index  which  counts  the  number  of  cycles.  It  is  incremented 
by  1  at  the  end  of  each  cycle. 

V(J)  =  specific  volume  for  Jth  increment  in  stress 
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EPX  ■  ,x 
EPY  -  Ey 
EXE 

EVE  '  */ 

EPXP  =  e  P 

A 

EPYP  =  e  P 
GAMME  =  Ye 
S(J)  =  entropy 

K  =  isothermal  bulk  modulus  =  K 

MU  =  y 

LAM  =  X 

YB  =  Y 

WPB  =  Wp 

TB  =  T  These  quantities  indicated  by  B  or  —  are  averages  of  values  at 
the  last  J  and  the  next  J.  Their  use  makes  the  integration 
scheme  approximately  second  order. 

DY  =  dY 

ISTOP,  IFLAG  are  control  indices  used  for  stopping  the  calculation  and  cal¬ 
culating  the  end  of  the  cycle. 

IY  is  an  index  used  to  guide  the  computation  through  the  elastic  branch  (IY=1) 
or  elastic-plastic  branch  (IY=2). 

JA  =  1  indicates  region  FOA  of  Fig.  4. 

=  2  indicates  region  AB  of  Fig.  4. 

=  3  indicates  region  BC  of  Fig.  4. 

=  4  indicates  region  CF  of  Fig.  4. 

OM  is  the  index  which  directs  the  computation  to  the  next  branch  of  the  cycle. 
IS  is  an  index  used  to  calculate  the  end  points  (A,  B,  C,  F,  0)  of  each  branch. 
GAM  =  r 
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The  program  contains  a  subroutine  for  use  with  a  CALCOMP  plotter.  It 
is  by-passed  in  this  listing,  but  it  can  be  engaged  by  removing  card  number 
0233,  "60  T0  99." 

Following  the  listing  is  sample  output  for  the  two  cases  shown  in 
Table  D-l  for  isentropic  and  adiabatic  compression. 
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PROGRAM  "CYCLE" 


FORTRAN  IV  G  LEVEL  20 


MAIN 


DATE  *  72332 


15/18/08 


0001 
0002 
0C03 
0004 
0005 
COO  6 
0007 
0008 
0009 
0010 


00 

co 


DIMENSION  SIGX(1500)»SIGY( 1500 1 ,EPX( 1500) , EPYI1500) ,EXEl 1500)  ROOl 

D I  ME NS  I  ON  EYE  Cl  500) ,EPXPI  1503) , EPY P ( 1 500) , WP ( 1500) ,  Y( 1500)  R002 

DIMENSION  P(i500),V(15C0),JH  1500)  R003 

DIMENSION  T(IOOO) 

DIMENSION  S(1000)»  GAMME11000) 

REAL  K,  KO,  MU,  LAM,  MUO,  MU1 ,  MU2 

REAL  J 1 , LENGX  ,  LENGY  R005 

INTEGER  PENP  P.005AAA 

DATA  XEP  X  / '  EPX '  /  ,  XEXE / '  EXE'/.XEYE/*  EYE'/  P.0C5AA 

DATA  XJ1/'  J  */,XSIGX/'SIGX'/,XSIGY/'SIGY'/,XEPXP/'EPXP'/  R005A 


C 

C  KO,K,LAM,MU,MUO  ARE  IN  KILOBARS 

C  B  IS  DIMENSIONLESS 

C  BETA  IS  IN  KILOBARS/DEGREE  C,  BEAT=CV* IGAM0+GAM1 )/V00 
C  C V  IS  IN  KB-CC/GRAM-DEGREE  C, 

C  VfVOO  ARE  IN  CC/G 
C  GAM, GAMO, GAM1  ARE  DIMENSIONLESS 
C  DELTA  IS  IN  KILOBARS/DEGREE  C 
C  DELTA=- ( AL PHA*KO*B-BET A ) 

C  ALPHA=  VOLUME  EXPANSION  COEFFICIENT, (DEG  C > 75E-4  FOR  AL 
C  K= ISOTHERMAL  BULK  MODULUS 
C  LAM  AND  MU  ARE  LAME*  CONSTANTS 

C  SGAM  IS  IN  IDEG  C)**-l,  SGAM= ( 4*BETA/3K0-BETA/ ( 3*K0*B) -ALPHA* ( BETA ) **2 
C  IY=T  IF  COMPRESSION  10R  EXPANSION)  IS  ELASTIC 
C  I Y=2  IF  COMPRESSION  IS  ELASTIC-PLASTIC 

C  IC=1  FOR  ISOTHERMAL  PROCESS 

C  IC=2  FOR  I SENTROPIC  PROCESS 

C  IC=3  FOR  ADIABATIC  PROCESS 

C  CONSTANTS  R006 

C  ALUMINUM  R007 

C 


0011 

0012 

0013 


IC  =  1 
N=  100 
MM=2 


R010 


RO 15 
RO 16 
R017 


R019 

RC20 

R02 1 
R022 

R023 
R024 
R025 
R026 
R027 
R02  8 
R029 
RC30 

R031 


R035 


R032 


0048 

0C49 

C050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 


0059 

0060 


0061 
0062 
0063 
0064 
0065 
0066 
0C67 
0068 
0069 
0070 
CO  7 1 
0072 
0073 
0074 
CC75 
0076 
0077 
0078 


K=KO 
MU=MUO 

LAM=K-I2./3.)*MU 
YB  =  YO 
WPB=WP( 1) 

TB=T  l  1 ) 

DY=0. 

VB=VOO 
I ST0P=1 
IFLAG=1 
400  CONTINUE 
C 

C  ALONG  OA  R036 

WRITE(6fl0l) 

10  l  FORMAT (4X»tJ*t8X»,SIGX*»llX*,SIGY,»13Xt ’P* , 14X , »Y» , 14X, 'VSllX, 

♦  •EPXSllX*  'GAMME  • ) 

C 

WR I TE l 6» 102 ) 

J=1  R037 

301  CONTINUE 
IY=l 

JA  =  l 
JM-2 
IS=1 

GO  TO  201 

302  CONTINUE 
I  Y  =  2 
JA=2 

JM  =  3 
IS  =  l 

GO  TO  201 

303  CONTINUE 
I  Y=  1 
JA=3 
JM«4 


00 

<T> 


0079 

DS  IGX=-DSIGX 

0080 

I  S=1 

0081 

GO  TO  201 

0082 

304 

CONTINUE 

0C83 

IY®2 

CC84 

JA-4 

CC85 

JM=l 

0086 

IS  =  1 

0087 

DS1GX=-0SIGX 

0088 

SIGXT=0.0 

0089 

201 

CONTINUE 

0090 

1 

SI  GXI  J+i)=SIGXI  J)«-DSIGX 

c 

IF(IFLAG.EQ.2.AND.SIGX(J«-1).LT*0.0)  GO  TO  2 
GO  TO  3 

2  CONTINUE 

DSI GX-— S I GXI J ) 

IFLAG=  1 
I S  TQP=2 
WRITE(6,110) 

110  FORMAT I2X, 'NEXT  VALUES  COMPLETE  THE  CYCLE') 
GO  TO  1 

3  CONTINUE 

C  COMPUTE  SIGYIJ+1) 

C 


0101 

DO  8  M= 1 » 3 

0102 

GO  TO  141,42), IY 

0103 

41 

OS  IGY=LAM*D$IGX/(LAM+2.*MU) 

0104 

GO  TO  43 

0105 

42 

DSIGY=DSIGX+DY 

0106 

IFIJA.EQ-4)  OSI GY=DSI GX— OY 

0107 

43 

CONTINUE 

01C8 

SIGYIJ  +  1)  =  SIGYI J  KDSIGY 

0109 

PIJ  +  1)=-ISIGXIJ*-1)+2.*SIGY<J  +  1  ))/3 

0110 

PB=(P ( J+l ) «-P( J) )/2. 

0111 

op*pi  j+d-pi  j) 

0112 

K=KO*B*PBMDELTA*SGAM*PB)*ITB-TOI 

0091 

C092 

0C93 

0094 

C095 

0096 

0097 

0C98 

0099 

0100 


R039 


R041 


R046 


*k  ft 


I 


0113 

DV*-VB*DP/K 

R048 

0114 

DEPX=-DP/K 

R047 

0115 

v(  j+n=vt  jh-dv 

R049 

0116 

VB= VI J ) +DV/2 • 

0117 

MU=MUO+AMU*( TB/VB-TO/VOO) 

0118 

r 

LAM=K-(2./3. )*MU 

1 

w 

c 

r 

COMPUTE  PLASTIC  WORK  ANO  Y 

0119 

V 

GO  TO  (5i»52)fIY 

0120 

51  DWP=0. 

0121 

DY=0. 

0122 

GO  TO  53 

0123 

52  CONTINUE 

0124 

IF( JA-EQ-4) YB=-YB 

0125 

I F I JA.E0.4IDY=-DY 

0126 

DWP=-I 2.*VB*YB/3.I*(DEPX+0Y/C2**MU)) 

0127 

DY=(A/WO)*DWP*EXPI-WPB/WO) 

0128 

53  CONTINUE 

0129 

WP(  JM  )=WPt  JH-DWP 

R086 

0130 

WPB  =  (  WPI  J+n*WP("J>  1/2, 

R087 

0131 

Y«  J*1I»Y4  J)«-DY 

R085 

0132 

r 

YB=(Yl  J«-1)»Y(  J))/2. 

R140 

c 

c 

COMPUTE  GAMMA-ELASTIC 

0133 

V 

DGAMME=(DSIGX-DSIGY)/(4.*MU1 

0134 

GAMMEt  J  +  l )=GAMME( J) +DGAMME 

0135 

GAMMEB^GAMME  I  JKDGAMMF/2. 

c 

C  COMPUTES  S  AND  T 


0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 


0151 


0152 

0153 

0154 

0155 

0156 


0157 

0158 

0159 

0160 

0161 

0162 


C 

GAM=GAM0*GAM1*< 1.-V8/VOOJ 
GO  TO  (31 *32 *33) » IC 

31  DT=0 

DS-GAM*CV*DV/VB-l6.*GAMMEB*AMU*DGAMME/3. 

GO  TO  34 

32  DS=0 

DT=-GAM*TB*DV/VB+16.*GAMMEB*TB*AMU*0GAMME/(3.*CV) 

GO  TO  34 

33  DS=DWP/TB 

DTXTB*DS/C V-GAM*TB*DV/  VB+ 16**3 AMMEB*T  8* AMU*OGAMME/  ( 3»*CV 1 

34  T( J+l ) =TI J ) *DT 
S( J+1)=S( J)+OS 
TB*TI  J  )*DT/2* 

102  FORMAT (i5X,'EPXP%12X,  •  EPXE » 8  10X,  •  EPYE*  9 13X,  ' S  •  , 15X ,  •  WP •  ,10X ,  •  T»  1 
8  CONTINUE 
C 

c 

GO  TO  (61 «  62 ) • I Y 
C 

C  COMPRESSION  IS  ELASTIC 
C 

61  DEPXP=0. 

0EPXE=D6PX 

DEPYP=0. 

DEPYE=Q. 

GO  TO  63 

C 

C  COMPRESSION  IS  ELASTIC-PLASTIC 
C 

62  CONTINUE 

DEPXP=(2./3«)*( DEPX*DY/(2«*MU) J 

OEPXE=DEPX-DEPXP  «  R093 

DEPYP=— DEPXP/2* 

OEPYE=— OEPYP 

63  CONTINUE 


\ 


r} 


0163 

EPXI  JM)=EPXI  J)*DEPX 

0164 

EPXPt  J*l) =EPXP( JJ+DEPXP 

0165 

EXElJ+i)=EXE(J)+DEPXE 

0166 

EPYPtJ+l)=EPYP(J)+DEPYP 

0167 

EYEIJf  l)  =  EYE(  JH-DEPYE 

0168 

J=J+1 

0169 

WRITEI6, 104) J,SIGX(J) , SIGYI J)  ,  P ( 4 ) * Y( J ) , V {  J ) , EPXC J ) ,GAMME( J) 

0170 

104 

FORMAT (2X, I 4 ,5F1 5. 8 ,2E 15.6 ) 

0171 

WRITEl6,105)EPXPIJ),EXE(J) ,  EYE (J),S(J),WP{J) ,T(J) 

0172 

105 

F0RMAT(i0X,5El5.6,F15.8> 

0173 

I F I J • GT.999 ) GO  TO  21 

0174 

IF(IS.EQ.2)DSIGX=SIGXB/Q 

0175 

IF (IS.EQ.2)  GO  TO  76 

0176 

GO  TO  171,72,73,74) ,JA 

0177 

71 

CONTINUE 

0178 

IFIISTOP.EQ.DGO  TO  5 

0179 

IF  (MCOUNT.GT.MM)GO  TO  21 

0180 

4 

CONTINUE 

0181 

J  =  l 

0182 

I  STOP  =  1 

0183 

IFLAG=1 

0184 

IC*IC+1 

0185 

IFI  IC.GT.3 )G0  TO  21 

0186 

DS  IGX=SIGXB/Q 

0187 

GO  TO  400 

0188 

5 

CONTINUE 

0189 

IF  (  (  SIGXI  Jl-SIGYI  jn.GT.~YUH  GOTO  201 

0190 

ZN=S!GX( J)-SIGY( J) 

0191 

ZO=SIGX{ J-l)-SIGY(J-l) 

0192 

OZ=ZN-ZO 

0193 

DSIGX=-(DSIGX/DZ)*(Y<  J-1)*Z0) / I l.-DY/DZ) 

0194 

J  =  J-1 

0195 

ISMS  +  l 

0196 

WRI TE (6, 106 ) 

0197 

106 

FORMATUOX,  'NEXT  VALUES  ARE  FOR  POINT  AM 

0198 

GO  TO  201 

0199 

72 

IF(SIGXU).GT.SIGXB)  GO  TO  201 

0200 

DSIGX*SIGXB-SIGX( J-l) 

R091 

R090 

R094 

R096 

R097 


0201 

J=  J-l 

0202 

IS*IS«-1 

0203 

WRITE(6,107) 

0204 

107 

FORMAT ( 10X»  *  NEXT  VALUES  FOR  POINT  Q®  1 

0205 

GO  TO  201 

0206 

73 

IFUSIGX1  Jl-SlGYt  J)  l.LT.Yl  Jl)  GO  TO  201 

0207 

ZN=SIGX< Jl-SIGYt J) 

0208 

ZO«SIGX(  J-ll-SIGY!  J-l) 

0209 

DZ*ZN-ZO 

0210 

DSIGX=  IDSIGX/DZ)*(Y(  J-ll-ZOl/I  l.-DY/DZ) 

0211 

J  =  J-1 

0212 

I S* I S+ 1 

0213 

WRITE(6,108) 

0214 

108 

FORMA TIIOX^'NEXT  VALUESFOR  POINT  C«) 

0215 

GO  TO  201 

0216 

74 

IFIEPX(J) • l T.0.0 )  GO  TO  201 

0217 

IF1EPX1  J-U.LT.O.O)  D£PX*EPX(  J  I 

0218 

DSIGXT=(LAM+2.*MUJ*DEPX 

0219 

SIGXT=SIGXT*DSIGXT 

0220 

IFISIGXT.LT. SIGXC J) 1  GO  TO  201 

0221 

DSIGX=(OSIGX/OSIGXn*I SIGX1 J-l  1-SI GXT+DSIGXT)/ 1 i.-DSIGX/DS IGXT 1 

0222 

J  =  J-1 

0223 

IS=IS*1 

0224 

WRITE! 6f 109) 

0225 

109 

F0RMAT(10X»'NEXT  VALUES  FOR  POINT  F® ) 

0226 

IF  L AG  =  2 

022  7 

MC  OUN  T=MC0UNT+1 

0228 

GO  TO  201 

0229 

76 

CONTINUE 

0230 

GO  TO  ( 301 » 302  » 303  #304 1 »  JM 

0231 

21 

CONTINUE 

0232 

55 

CONTINUE 

0233 

GO  TO  99 

0234 

56 

CONTINUE 

0235 

DO  20  K2=1#J 

R160 


R162 


R163 


0236 

20  J1(K2)=K2 

R164 

0237 

ZER0= 1*0 

R165 

0238 

LENGX=4.5 

R 166 

0239 

LENGY=LENGX 

R 167 

0240 

NCHAR  X=4 

R 168 

0241 

NCHARY=NCHARX 

R 169 

0242 

PL  END  =  7  *0 

R170 

0243 

PE  NP=-3 

R 1 71 

0244 

OR  IGI N=0.0 

Rl  72 

0245 

FACT=i.O 

P  173 

0246 

CALL  PLOTITl J1 ,S I GX ,LENGX 

,LENGY,NCHARX,NCHARY,XJ1 

,XSIGX,PLEND,ZER0R174 

C,J,PENP, ORIGIN, FACT) 

R 1 75 

0247 

CALL  PLOTITt J1  ,S IGY , LENGX 

»  LENGY, NCHARX, NCHAR Y, XJ1 

,XSIGY,PL£ND,ZER0R176 

C, J» PE NP» ORIGIN, FACT) 

R177 

024  8 

CALL  PLOTI T{ Jl, EPXP, LENGX 

,LENSY, NCHARX, NCHARY, XJl 

, X EPXP, PLEND,ZEROR 178 

C, J,PENP,ORIGIN»FACT) 

R  179 

0249 

CALL  PLOTITI J1,EPX, LENGX, 

LENGY, NCHARX, NCHARY,XJi, 

XEPX,PLEND,ZERO, JR180 

C»PENP, ORIGIN, FACT) 

R 1  81 

0250 

CALL  PLOT  IT  1 Jl , EXE  * LENGX, 

LENGY, NCHARX, NCHARY, XJl, 

XEXE,PLEND,ZER0,JR182 

C,PENP, ORIGIN, FACT) 

Rl  83 

0251 

PE  NP=999 

R184 

0252 

PLEND= 10*0 

R185 

0253 

CALL  PLOTITI Jl, EYE, LENGX, 

LENGY , NCHARX , NCHARY, XJl, 

XEYE, PL  END, ZERO, JR  186 

C,PENP,0RIGIN,FACTI 

Rl  87 

0254 

99  STOP 

Rl  88 

0255 

END 

R189 
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SUBROUTINE  PLOT! T f X, Y , L ENGX,LE NGV, NCHARX,  NCHAR Y,TITLEX, TITLEY* PLENS100 


CD,ZERQ»NPTS*PENP»QRIGN*FCTR1  SlOl 

DIMENSION  IBUF(1500),XU500J,VC  1500)  S102 

REAL*8  TITLEX, TITLEY  S103 

REAL  LENGX#LENGY  S104 

INTEGER  PENP  S105 

DATA  ENTER/1./  S106 


INSTRUCTIONS  FOR  USING  THE  PLOT  SUBROUTINE; 

THE  TITLES  FDR  THE  PLOTS  MUST  BE  LITERAL  CONSTANTS  NO 
LONGER  THAN  EIGHT  CHARACTERS  IN  LENGTH  AND  MUST  BE  DECLARED 
REAL*8. 

FOR  THE  CASE  OF  TOTALLY  POSITIVE  DATA  THE  SCALE  ROUTINE 
WILL  RETURN  A  NON-ZERO  ORIGIN  FOR  ONE  OR  BOTH  AXES.  IF  A 
ZERO  ORIGIN  IS  DESIRED  FOR  BOTH  AXES  THE  VARIABLE  ORIGN  SHOULD 
BE  SET  TO  ZERO  AND  FOR  NON-ZERO  ORIGINS  ORIGN  SHOULD  BE  SET  TO 
ONE. 

CONTINUE 

IN  CALLING  THE  SUBROUITNE  THE  VALUE  ASSIGNED  TO  THE  VAR¬ 
IABLE  ZERO  SHOULD  BE  NO  SMALLER  THAN  .5  SINCE  ANYTHING  SMALLER 
WOULD  NOT  ALLOW  ROOM  TO  PRINT  THE  TITLE  FOR  THE  X-AXIS. 

CARE  MUST  BE  TAKEN  IN  CHOOSING  THE  VALUES  FOR-  THE  VAR¬ 
IABLE  LENGY.  THE  SUM  OF  LENGY  AND  ZERO  SHOULD  NOT  EXCEED  10. 
THE  WIDTH  OF  THE  PLOTTING  PAPER  IS  ONLY  11  INCHES  AND  THE 
VALUE  FOR  ZERO  TELLS  THE  PLOTTER  HOW  FAR  ABOVE  THE  BOTTOM  OF 
THE  PAPER  TO  BEGIN  DRAWING  THE  Y-AXIS,  THEREFORE  THE  AXIS  MUST 
BE  SHORT  ENOUGH  TO  FIT  ON  THE  PAPER  WITHOUT  HITTING  THE  LIMIT 
SWITCH  ON  THE  TOP. 

THE  VARIABLE  FCTR  IS'  USED  TO  CHANGE  THE  SIZE  OF  THE  PLOTS 
WITHOUT  CHANGING  THE  LENGTH  OF  THE  AXES.  IT  PRODUCES  A  PLOT 
WHICH  IS  A  SCALEO  REPRESENTATION  OF  THE  NORMAL  PLOT.  IF  FCTR 
IS  NOT  EQUAL  TO  ONE  CARE  MUST  BE  TAKEN  TO  INSURE  THE  NEW  PLOT 
WILL  NOT  BE  TOO  LARGE  TO  FIT  ON  THE  PAGE  OR  TOO  CLOSE  TO  THE 


o 


i*  fJ 


CC08 


>x> 

CO 


0C09 


* 


C  BOTTOM  LIMIT  SWITCH.  A  PLOT  MADE  VERY  SMALL  BY  THE  USE  OF  THE 

C  FCTR  PARAMETER  MAY  NOT  BE  READABLE,  THEREFORE  THE  USE  OF  FCTR 

C  SHOULD  BE  LIMITED. 

C  THE  VARIABLE  PLEND  IS  THE  NUMBER  OF  INCHES  FROM  THE 

C  ORIGIN  THE  PEN.  WILL  MOVE  IN  THE  X  DIRECTION  AT  THE  END  OF  THE 

C  PLOT.  THEREFORE  THE  VALUE  ASSIGNED  TO  PLEND  SHOULD  BE  THE 

C  LENGTH  OF  THE  X-AXIS  PLUS  THE  NUMBER  OF  INCHES  PAST  IT  THE 

C  NEW  PLOT  IS  TO  BEGIN.  THE  VALUE  ASSIGNED  TO  PLEND  IN  THE 

C  FINAL  CALL  SHOULD  BE  AT  LEAST  THE  VARIABLE  LENGX  +  5  TO  ALLOW 

C  THE  PAPER  TO  BE  TORN  OFF  THE  PLOTTER  WITHOUT  AFFECTING  THF 

C  PLOT. 

CONTINUE 

C  THE  VARIABLE  PENP  SHOULD  BE  ASSIGNED  THE  VALUE  -3  FOR 

C  ALL  CALLS  TO  THE  SUBROUTINE  EXCEPT  THE  FINAL  CALL.  FOR  THE 

C  FINAL  CALL  THE  VALUE  999  SHOULD  BE  GIVEN  TO  PENP.  THESE 

C  VALUES  MERELY  RAISE  THE  PEN  WHILE  MOVING  TO  THE  NEW  ORIGIN, 

C  OR  FOR  THE  CASE  WHEN  PENP=999,  WHILE  MOVING  TO  THE  END  OF  THE 

C  PAPER. 

C 

CONTINUE 

C 

C  Xs  ARRAY  CONTAINING  X  COORDINATES 

C  Y=  ARRAY  CONTAINING  Y  COORDINATES 

C  LENGX=  LENGTH  OF  X-AXIS 

C  LENGY=  LENGTH  OF  Y-AXIS 

C  NCHARX=  NUMBER  OF  CHARACTERS  IN  TITLE  OF  X-AXIS 

C  NCHARY  =  NUMBER  OF  CHARACTERS  IN  TITLE  OF  Y-AXIS 

C  T I T  LE  X  =  TITLE  TO  BE  GIVEN  X-AXIS 

C  TITLEY=  TITLE  TO  BE  GIVEN  Y-AXIS 

C  PLEND*  DISTANCE  PEN  WILL  MOVE  AT  END  OF  PLOT 

C  ZERO-  DISTANCE  ABOVE  BOTTOM  OF  PAGE  PEN  SHOULD  BE  MOVED  TO 

C  BEGIN  Y-AXIS 


0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
C02 1 
0C22 
002  3 
0C24 
002  5 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


NPTS*  NUMBER  OF  POINTS  ♦  2  IN  ARRAYS  TO  BE  PLOTTED  (  NUMBER 
OF  X  POINTS  SHOULD  EQUAL  NUMBER  OF  Y  POINTS.  THE  TWO 
EXTRA  PLACES  IN  THE  ARRAYS  ARE  USED  IN  THE  PLOTTING 
ROUTINE  FOR  SCALING) 

PENP=  POSITION  OF  PEN  AT  END  OF  PLOT  (NORMALLY  EITHER  -3  OR 
999) 

OR  I GN=  0.0  FOR  ZERO  ORIGINS  AND  1.0  FOR  NON-ZERO  ORIGINS 

FCTR=  RATIO  OF  NEW  PLOT  SIZE  TO  THE  NORMAL  PLOT  SIZE.  A  VALUE 
OF  1.0  GIVES  PLOTS  OF  NORMAL  SIZE.  ANYTHING  GREATER 
THAN  ONE  INCREASES  THE  SIZE  OF  THE  PLOT  WHILE  ANYTHING 
LESS  THEN  DECREASES  THE  PLOT  SIZE 


IF (ENTER. NE.l)  GO  TO  1 
CALL  PLOTS! I BUF, 1000) , 
l  ENTER=2.0 

CALL  PLOTIO. 0,-11. 0,-3) 

CALL  FACTOR(FCTR) 

CALL  PLOT  10. 0, ZERO, -3) 

CALL  SCALE(X,LENGX,NPTS,1) 

CALL  SCALE ( Y, L ENG Y » NPTS ,1) 

XPAGE=0.0 

YPAGE=0.0 

I  FI YtNPTS*l)  )  10,20,30 

20  IF ( X ( NPTS* 1 ) )  21,22,23 

22  CALL  AXIS (0.0, YP AGE , T I TLEX ,-NC HARX, LENGX, 0.0, X (NPTS  +  1 ) , X (NPTS+2  ) 
CALL  AXIS! XPAGE, 0.0, TI TLEY , NCH ARY, LENGY ,90 .0,  Y ( NPTS*1 ) 9 YI NPTS+2 ) 
CALL  LINE!  X , Y, NPTS, 1 , 0 ,0) 

CALL  PL0T(PLEND,-11.0, PENP) 

GO  TO  99 

21  XPAGE=-X(  NPTS+D/X  (NPTS+2) 

GO  TO  22 

23  IF (ORIGN.NE.O.O)  GO  TO  22 
XMAX=X( NPTS*1 )*LENGX*X(NPTS*2) 

X( NPTS+1)=0.0 

X( NPT  S+2) “XMAX/LENGX 
GO  TO  22 
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S1C9 

5110 
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5111 
SI  12 
SI  13 
S113A 
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)  S116 
)  SI  17 
S118 
SI  19 
$120 
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S123 

5127 
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0034 

0035 

0036 

0037 

0038 

0039 

0C40 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 


10  YPAGE=— YCNPTS+1)/Y{NPTS+2J 
I F l X { NPTS+l ) )  11,12,13  * 

12  GO  TO  22 

11  XPAGE=-XI NPTS+ll/X(NPTS+2) 

GO  TO  22 

13  IF(ORIGN.NE.O.O)  GO  TO  12 
XMAX=X(NPTS+l)+LENGX*X(NPTS+2) 
XI NPT  S  +  i } -0»0 

X ( NPT  $+2)=XMAX/LENGX 
GO  TO  22 

30  IFIXINPTS+i)  »  31,32,33 
32  IF(OR IGN.NE.C)  GO  TO  22 
34  YMAX= Y< NPTS+l )+LENGY*Y(NPTS+2) 
Yl  NPTS  +  l )  =  0. 0 
Y( NPTS+2) =YMAX/LENGY 


5132 

5 133 

5135 

5136 

5137 

5138 

5139 

5140 

5 14 1 

5142 
SI  58 
SI  59 

5160 

5161 

5162 


0049  GO  TO  22  S163 
0050  31  XPAGE--XC NPTS+l ) /X { NPTS+2 )  S164 
0051  GO  TO  32  S166 
0052  33  IFtORIGN.NE.O.O)  GO  TO  22  S170 
0053  XM AX= XI NPTS+l )+LENGX*X I NPTS+2)  S171 
0054  X ( NPTS+l ) =0«0  S172 
0C55  X(NPTS+2)*XMAX/LENGX  S173 
0056  GO  TO  34  S174 
0057  99  RETURN  S175 
0058  .  END  S176 


J  SIGX  SIGY  P  Y  V  EPX  GAMME 

EPXP  EPXE  EPYE  S  WP  T 

2  -0.40599990  -0.20910841  0.27473891  10.00000000  0.36986703  -0. 359291 E-03  -0. 179646E-03 

0.0  -0.359291 E-03  0.0  -0.603608E-05  0.0  300.00000000 

3  -0.81199980  -0.41821682  0.54947758  10.00000000  0.36973411  -0 . 71 8581 E-03  -0.359291E-03 

0.0  -0. 718581E-03  0.0  -0. 120722E-04  0.0  300.00000000 


51-20.29986572  -10.45541763  13.73689747  10.00000000  0.36341119  -0. 179644E-01  -0.898225E-02 

0.0  -0. 179644E-01  0.0  -0. 301802E-03  0.0  300.00000000 

52  -20.70585632  -10.66452599  14.01163197  10.00000000  0.36328059  -0. 183237E-01  -0.916189E-02 

0.0  -0. 183237E-01  0.0  -0.307838E-03  0.0  300.00000000 

NEXT  VALUES  ARE  FOR  POINT  A 

52  -20.62063599  -10.62062836  13.95396423  10.00000000  0.36330801  -0. 182483E-01  -0. 912418E-02 

0.0  -0. 182483E-01  0.0  -0.306570E-03  0.0  300.00000000 

£  53-21.02662659  -11.02662754  14.35995960  10.00000000  0.36311513  -0. 187792E-01  -0. 912418E-02 

-0.353961 E-03  -0.184252E-01  -0. 1 76981 E-03  -0.315490E-03  0.128563E-02  300.00000000 


101  -40.51417542  -30.51428223  33.84757996  10.00000000  0.35397691  -0.442640E-01  -0.912418E-02 

-0.173439E-01  -0.269202E-01  -0.867194E-02  -0.743629E-03  0.621992E-01  300.00000000 

102  -40.92016602  -30.92027283  34.25357056  10.00000000  0.35378897  -0.447950E-01  -0. 912418E-02 

-0.176978E-01  -0. 270972E-01  -0.884892E-02  -0. 752549E-03  0.634518E-01  300.00000000 

NEXT  VALUES  FOR  POINT  B 

102  -40.59999084  -30.60009766  33.93339539  10.00000000  0.35393715  -0.443763E-01  -0. 912418E-02 

-0. 174187E-01  -0.269576E-01  -0.870935E-02  -0.745514E-03  0.624640E-01  300.00000000 

103  -40.19398499  -30.39097595  33.65864563  10.00000000  0.35406429  -0.440170E-01  -0.894453E-02 

-0-1741 87E-01  -0.265983E-01  -0.870935E-02  -0. 739478E-03  0.624640E-01  300.00000000 
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p 

EPYE 


Y 


EPX 

T 


GAMME 
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204 


SIGX 

EPXP 


SIGY 

EPXE 


V 

WP 


0.40636444  -9.47920322  6.18401051  10.00000000  0.36701411  -0.808666E-02  0.902007E-02 

-0.1741 87E-01  0.933193E-02  -0.870935E-02  -0.135848E-03  0.624640E-01  300.00000000 

0.81236434  -9.27009392  5.90926552  10.00000000  0.36714596  -0.772735E-02  0.919972E-02 

-0 . 1741 87E-01  0.969123E-02  -0.870935E-02  -0.12981 2E-03  0.624640E-01  300.00000000 


NEXT  VALUES  FOR  POINT  C 


204  0.64233136  -9.35766888  6.02433205  10.00000000  0.36709070  -0.787783E-02  0.912448E-02 

-0 . 1741 87E-01  0. 954075E-02  -0.870935E-02  -0. 1 32340E-03  0.624640E-01  300.00000000 

205  1.04833126  -8.95166874  5.61833191  10.00000000  0.36728561  -0. 734688E-02  0.912448E-02 

-0. 170647E-01  0.971773E-02  -0.853236E-02  -0. 123420E-03  0.637637E-01  300.00000000 


253  20.53619385  10.53631401  -13.86960316  10.00000000  0.37676501  0.181 385E-01  0.912448E-02 

-0.743989E-04  0.1821 28E-01  -0. 371 640E-04  0.304734E-03  0.126967E  00  300.00000000 

254  20.94218445  10.94231319  -14.27559853  10.00000000  0.37696511  0.186694E-01  0.912448E-02 

0.279562E-03  0.183898E-01  0.139817E-03  0.313654E-03  0.128301E  00  300.00000000 

NEXT  VALUES  FOR  POINT  F 

254  20.61920166  10.61933422  -13.95262146  10.00000000  0.37680590  0.182470E-01  0.912448E-02 

-0.202062E-05  0.182490E-01  -0.974804E-06  0.306558E-03  0.127240E  00  300.00000000 

255  20.21319580  10.41022491  -13.67788124  10.00000000  0.37667054  0.178877E-01  0.894483E-02 

-0. 202062E-05  0.178897E-01  -0.974804E-06  0.300521E-03  0.127240E  00  300.00000000 


303  0.72513318  0.37298089  -0.49036467  10.00000000  0.37022889  0.641669E-03  0.321765E-03 

-0. 202062E-05  0.643598E-03  -0.974804E-06  0.107875E-04  0.127240E  00  300.00000000 

304  0.31913328  0.16387248  -0.21562606  10.00000000  0.37009585  0.282378E-03  0. 1 421 1 9E-03 

-0.202062E-05  0.284308E-03  -0.974804E-06  0.475146E-05  0.127240E  00  300.00000000 

NEXT  VALUES  COMPLETE  THE  CYCLE 

305  0.0  -0.00049567  0.00033045  10.00000000  0.36999130  -0. 393484E-07  0.9101 79E-06 

-0.202062E-05  0.189035E-05  -0. 974804E-06  0. 684031 E-08  0.127240E  00  300.00000000 


APPENDIX  E 


Gibbs  Relation 

Consider  an  element  of  mass  within  which  may  operate  various  irreversible 
processes.  The  element  is  supposed  to  be  in  a  stable  or  equilibrium  state;  if  it 
is  left  undisturbed  it  remains  unchanged.  We  also  suppose  that  temperature  T, 
specific  entropy  S,  and  specific  internal  energy  E  can  be  assigned  to  it.  External 
forces  are  applied  at  a  particular  instant  and  work  is  done  on  the  element  in  the 
amount  dW.  Heat  of  amount  dQ  is  also  added.  These  external  influences  set  in 
motion  irreversible  processes,  some  of  which  operate  immediately  and  some  of  which 
may  operate  over  an  extended  time  interval.  After  a  long  time,  during  which  the 
element  is  insulated  from  further  heat  exchange,  the  element  is  assumed  to  come  to 
a  new  equilibrium  state  which  differs  from  the  first  by  internal  energy  dE, 
temperature  dT  and  entropy  dS.  We  suppose  that  the  only  net  effect  of  adding  dQ 
and  dW  is  to  increase  the  internal  energy: 

dE  =  d Q  +  dW  (El) 

We  further  suppose  that  some  method  exists  for  assessing  the  internal  changes  in 
the  sample,  i.e.  for  partitioning  the  change  dE  into  changes  in  thermal  energy, 

TdS,  and  coherent  energy,  say  dWg.  For  the  purposes  of  Section  IV  this  is 
accomplished  by  writing  the  total  work  increment  as  the  sum  of  an  elastic  and 
plastic  term  and  assuming  that  the  elastic  work  is  reversible: 

dW  =  dWe  +  dWp  (E2) 

dE  =  TdS  +  dWe  (E3) 

Combining  Eqs.  (El)  -  (E3)  yields  the  Gibbs  relation 

TdS  =  dWp  +  dQ  (E4) 

Properties  of  the  material  and  properties  postulated  to  this  point  imply 
that  there  also  exist  other  thermo-dynamic  functions,  e.g.  the  Helmholtz  free  energy, 


99 


A  =  E  -  TS 

dA  =  -SdT  +  dWQ 
e 

=  -SdT  +  dW  -  dW  (E5) 

P 


100 


APPENDIX  F 

Heat  Absorbed  in  Isothermal  Compression 
=  internal  energy  at  A 

Eg  =  internal  energy  at  B 

=  EA  +  dQ  +  dW 

dQ  =  heat  absorbed  by  sample  in  going  ' 

€x 

from  A  to  B  isothermal ly,  <0. 

dW  =  work  performed  on  the  sample  in  taking  it  from  A  to  B  isothermally,  >0. 
E  =  E(T ,  ex,  Ye) 

Ed  =  E.  +  ( 3E/3e  )  de  +  (3E/3Ye)  dYe 

B  T,ye  X  T,ex 

From  Eq.  (47) 

dE  =  TdS  -  Vpde  +  8VxdYe/3 

X 

dA  =  -SdT  -  Vpde  +  8VxdYe/3 

X 

Since  de  =  dV/V,  these  can  also  be  written 
✓\ 

dE  =  TdS  -  pdV  +  8VxdYe/3 
dA  =  -SdT  -  pdV  +  8VxdYe/3 

The  elastic  work  consists  of  two  terms:  one  arising  from  density  change  and 
one  from  elastic  shear.  Then 

Er  =  E»  +  (3E/3V)  p  dV  +  (9E/3Ye)  dYe 
B  A  T ,Ye  T,V 

(3E/3V)  _  =  (3E/3V)  p  +  (3E/3S)  (as/a V) 

T,Y6  S,Ye  V,Ye  T  ,Ye 

=  -p  +  T  (3S/3V) 

T,YS 
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S  =  -(3A/3T)  ,  p  =  -(3A/3V)  o 

V,Ye  T,Ye 

OS/3  V)  =  -(32A/3V3T)  =  (3P/3T) 

Ye  V,Ye 

Similarly, 

*  V 

(3E/3Y®)  .  =  (9E/3y6)  +  (3E/3S)  _  ( 9S/3Ye) 

T,v  S,V  V,Ye  T,V 


=  8Vt/3  +  T(3S/3Ye) 

T,V 


(3S/3y6)  =  -(32A/3T3Ye)w  =  -(8/3)( 3Vt/3T) 

T,V  V  V,Ye 

=  -(8V/3)  (3t/3T) 

V,Ye 

Then 

dQ  =  EB  "  EA  '  dW  =  £~p  +  t(3POT)  e  ]  dV 

V,y  • 


+  [8Vt/3  -  (8VT/3)  (3x/3T)  J  dYe 

V,Ye 


+  pdV  -  8VTdYe/3  -  8VTdYP/3 


dQ  =  T(3p/3T)  dV  -  (8VT/3)  (3t/3T)  dY6  -  8VrdYP/3 

V, y  V,Ye 


where  8VxdYP/3  =  dWp. 
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APPENDIX  G 


Evaluation  of  (9av/3T)  e  r  e  and  (9x/3T) 

x  ex  »  ey  Ye 

From  Eq.  (27) 

dax  =  (a  +  2 y)  de^Q  +  2a  deye 

where  A  and  y  are  functions  of  temperature  and  strain. 

Let  =  a(T2»e) 

A-|  =  A (T-j  , e ) 

y2  =  y(T2»e) 
y-j  =  y(T-|,e) 

0  0 

In  the  neighborhood  of  a  fixed  pair  of  elastic  strains,  Uxol»  eyol)>  °x  can  be 

.  e  6 

represented  as  a  linear  function  of  ox>  e  : 


a 


x 


=  °x^Tl’  £xol*  £yol ^ 

+  ( A-j  +  2y-| )  (exe  -  e®ol)  +  2A-|(eye  -  eyol ) 


Similarly 


°x(T2’  £x6’  ey6)  = 


xu29  xo2*  yo2 


»  ewn'?)  +  (^9  +  2y?)  (ex  -  Exo2) 


+  2A2  (ey  -  eyo2). 


Choose  ex02  and  £yo2  so  that  ax(T2’  ex02»  eyo2)  =  °x(Tl*  £x01*  eyol 


e  \ 
e..„n  ) . 


Then 


where 


6ox  =  (A2  +  2y2)  (e/  -  Exo2)  +  2A2(ey6  “  £yo2} 

-  (A1  +  2p 1 )  (ex  -  Gxoi)  ■  2Mey  ‘  £yol ^ 
6ax  —  ax(T2»  Ex  *  Gy  )  **  ^x  ( l"i  *  £x  9  Gy  ^  * 
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We  have  approximately 


*2  =  X-,  +  (3X/3T)6T,  v2  =  +  (3y/3T)6T 

exo2  ‘  exol  +  \  ST 
eyo2  =  £yol  +  “y  ST 


6aX  =  "aX^l  +  2 W])  6T  -  2xi  ay  6T 

+  6T(ex6  “  4l)  +  2y)  /  3T 

+  26T  (Eye  -  e*ol)  (3X/3T) 

(3ox/3T)  e  e  =  1  in)  (fio  /6T)  =  -  a  fx  +  2y)  “2X-i  a 

x  y  e  e  x  1  y 

ex  ^xol 


If  =  ay  =  a/3,  where  a  is  the  volume  expansion  coefficient,  then 

(3a  /3T)  e  e  =  -a(x  +  2y/3)  =  -  Ka 
X  ex  ,ey 

=  -(3p/3T)v 
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Let  dr  =  2ydye 

Let  t(T-|  ,y6)  =  t(T1  ,  yQle)  +  2y (T-j )  (ye  -  y^6) 

T(T2»ye)  =  t(T2>  yo2e)  +  2y  (T2 ) (ye  -  yQ2e) 

Let  t(T2,  yQ2e)  =  t(T-|,  y^6) 

Then 

t(T2»Y^)  -  t(T-|  >ye)  =  6t  =  2y(T2)  (  y6  -  Yo2e)  ~  2y (T-j )  (y6  -  y^6) 
Let  y(T2)  =  y(T-j)  +  (dT/dy)  (T2  "  T-j)* 

Then 

6t  =  -2y(T1 )  (yo2e  -  Yo1e)  +  2(T2  -  T, )  (dT/dy)  (Ye  -  y^). 

Let  6T  =  T2  -  T1  ->  0,  yo2e  -  y^6,  y6  -  Y^® 

Then  p 

6t/6T  ■+•  (  3t/  3T)  e  =  -2y(T1)(3ye/8T)T 

Y 

ye  =  const,  implies  dye  =  d(exe  -  de  e)/2  =  0. 
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APPENDIX  H 


Multiplication  Time  for.  Dislocations 

Suppose  that  each  dislocation  runs  freely  between  imperfections  without 
multiplying,  reaches  an  imperfection,  hesitates,  multiplies  and  goes  on  to  the 
next.  If  the  distance  between  imperfections  is  a  and  shear  velocity  in  the  crystal 
is  v$,  the  multiplication  time  is  the  order  of 

t0  -  »/»s  +  th 

where  t^  =  hesitation  time  at  the  imperfection. 

Let  t^  be  determined  by  a  rate  process.  A  dislocation  is  trapped  at  a 
site  and  must  surmount  an  energy  barrier,  U,  to  escape.  It  makes  n  escape 
tries  per  second  and  the  probability  that  any  one  of  them  will  succeed  is  assumed 
to  be  exp(-U/kT).  Let  P(M  +1)  =  probability  that  escape  occurs  on  the  (M  +  1; 

trial.  Then 

P(M  +  1 )  =  [1-  exp(-U/kT)]M  exp(-U/kT) 

=  (probability  that  first  M  tries  fail)  x 
(probability  that  (M  +  l)*-*1  succeeds). 

A1so  P(M  +  J)  =  P(M  +  1)  [1  -  exp  (-  U/kT)]0'1. 

The  probability  that  escape  occurs  on  the  (M  +  l)^*1,  (M  +  2)^ . .  or  (M  +  J)^ 

try  i  s 

J  i 

l  P(M  +  r)  =  P(M  +  1)  (1  -  a  )/ (1  -  a) 
r=l 

where  a  =  1  -  exp  (-U/kT). 

Since  there  are  n  tries  per  second,  M  =  nt.  Let  J  =  ndt, 
then 

1  -  ak  =  -  ndt  Jin  [1  -  exp  (-U/kT)]  +  0(n3t2) 
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Neglecting  the  higher  order  quantities,  we  have  the  result  that 

P(t)dt  =  -[1  -  exp(-U/kT)]nt  £n[l  -  exp(-U/kT)]  ndt 

where  P(t)dt  =  probability  that  escape  occurs  during  the  interval  (t,t  +  dt). 

It  is  readily  verified  that 

oo 

J0  P(t)  dt  =  1 

Then  the  mean  time  for  escape  is 

th  =  /“tP(t)dt 

=  -(1/n)  /  £n [1  -  exp(-U/kT)] 
a  exp(U/kT)/n  if  exp(-U/kT)  <<  1. 

The  multiplication  time  is  then 

tQ  =  £/vs  -  (1/n)  /  an[l  -  exp  (-U/kT)]  (HI) 

for  exp  (-U/kT)  <  1.  The  effect  of  applied  stress  is  to  diminish  U,  e.g. 

U  =  U0  (1  -  o/a0) 

where  aQ  is  a  critical  or  threshold  stress.  For  c  >  aQ,  Eq.  (HI)  is  replaced  by 

‘o  ■  l/s- 
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APPENDIX  I 

Functional  dependence  of  p  and  y 


32p/3T3Ye  =  -(8/3)  3(V3t/3T)/3V  (86a) 

32p/3N3T  =  -32r/3V3T  (86b) 

32r/3T3ye  =  (8/3)  3(V3t/3T)/3N  (86c) 

Let  t  =  2yye  (II ) 

3y/3ye  =  0  (12) 

Let  r  =  b2yV/2  (13) 

Eq.  (86a),  with  (II)  becomes 


32p/3T3ye  =  -(16ye/3)  3(V3y/3T)/3V 


=  - ( 1 6y 6/ 3 )  32(yV)/3T3V  (14) 

3p/3ye  =  0  (15) 

Then  Eq.  (14)  becomes 

32(pV)/3T3V  =  0  (16) 

or  3(uV)/3T  =  f-j  (T ,ye*N)  (17) 

In  view  of  Eq.  (12),  (17)  becomes 

3  (yV)/3T  =  f-j  (T , N) .  (18) 

From  Eqs.  (86b)  and  (13), 

32p/3N3T  =  - ( b2/ 2 )  32pV/3V3T  (19) 
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Equations  (19)  and  (16)  combine  to  yield 

32p/3N3T  =  0.  (no) 

Equation  (110)  implies  that 

3P/3T  =  f2(T,Ye.V)  (Ill) 

Since  p  is  already  assumed  to  be  independent  of  ye,  Eq.  (Ill)  becomes 

3P/3T  =  f2(TtV).  (112) 

Equations  (86c),  (II)  and  (13)  combine  to  give 

(b2/2)  32(yV)/3T3Ye  =  (16Ye/3)  32(yV)/3T3N  (113) 

Bv  Eq.  (12)  y  is  independent  of  ye>  so  the  left  hand  side  of  Eq.  (113)  is 
zero.  Then 

32(yV)/3T3N  =  0  (114) 


SO, 

3(yV)/3T  =  f3(T,V,ye)  (115) 

Comparing  Eqs.  (18)  and  (115)  shows  that 

3(yV)/3T  =  f i ( T )  (116) 

Integration  of  Eqs.  (112)  and  (116)  yields 

P  =  f(T,V)  +  g(V,N)  (117) 

yV  =  F(T)  +  G(V,N)  (118) 

where  the  assumptions  that  p  and  y  are  independent  of  y6  have  been  carried 
forward  into  Eqs.  (117)  and  (118). 

no 


APPENDIX  J 


Temperature  Dependence  of  Shear  Modulus  of  A1 

Equation  (75)  is  obtained  from  Eqs.  (70)  under  the  assumptions  that 

pressure  and  specific  heat,  C^,  are  both  independent  of  ye,  and  that  x  is 

proportional  to  ye: 

y  =  y0  +  b(T/V  -  T0/V0)  +  ,<j>(V)  (Jl) 

Such  a  linear  dependence  is  quite  reasonable  for  small  deviations  from  TQi 
but  must  clearly  be  inadequate  as  T  approaches  the  melting  temperature,  T  . 
If  one  assumes  that  the  reduction  of  y  with  temperature  is  an  activated  pro¬ 
cess  and  that  y  vanishes  at  T  =  Tm,  he  obtains  an  equation  of  the  form 

y  =  c  exp[-A/(Tm  -  T)]  (02) 

This  function  varies  slowly  with  T  until  it  approaches  Tm,  then  decreases 
precipitously  in  value.  To  give  it  the  linear  form  of  Eq.  (Jl),  for  small 
T  -  T  ,  let  y  be  given  by  the  expression 

y  =  c  exp[-A/(Tm  -  T)]  [1  -  b(T  -  TQ) ]  (03) 

Measurements  made  by  P.  M.  Sutton  (1953)  have  been  used  to  evaluate  the  con¬ 
stants  in  Eq.  (J3)  with  TQ  =  273. 2°K.  The  result  is 

y(megabars)  =  .2812  exp  [-7.8059/(933.2  -  T)] 

[1  -  4.79  x  10'4  (T  -  273.2)]  (J4) 

Values  calculated  from  Eq.  (J4)  are  compared  with  Sutton's  measured  values  in 
the  following  table: 
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T  ■  = 

273.2 

373.2 

473.2 

573.2 

673.2 

773.2  873.2  903.2 

933.2 

pcalc  = 

.2812 

.2672 

.2530 

.2384 

.2232 

.2061  .1780  .153 

0 

pmeas  = 

.2812 

.2670 

.2530 

.2385 

.2232 

.2063 

These  values  are  plotted  in  Fig.  Jl. 

It  would  improve  the  form  of  (J3)  to  write  it  as 

y  =  c  exp[-A/(Tm  -  T)]  [1  -  B'(T/V  -  T0/VQ)]  (05) 

In  that  case 

b'  -  1.8  x  10"4  cc/g°K  . 

However,  one  can  argue  that  the  parameter  A  should  also  be  volume  dependent, 
so  the  refinement  of  (J5)  hardly  seems  worthwhile  at  present. 

A  different  form  of  the  temperature  dependence  of  elastic  constants 
has  been  given  by  Varshni  (1970).  It  describes  the  slow  variation  satis¬ 
factorily  but  does  not  vanish  at  T  =  T  .  For  the  present  purpose  the 
behavior  near  T  =  Tm  is  important  because  it  relates  to  the  annealing 
properties  of  dislocations. 
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Figure  J1  C44  of  A1  according  to  the  formula 

C44  =  .2312  exp[-7. 3059/(933. 2-T] 

[1  -  4.79  x  10‘4  (T-273.2)] 

This  corresponds  closely  to  data  given 
by  P.  M.  Sutton  (1953). 
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